Dynamical spin-spin correlation functions in the Kondo model out of equilibrium 
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We calculate the dynamical spin-spin correlation functions of a Kondo dot coupled to two nonin- 
teracting leads held at different chemical potentials. To this end we generalize a recently developed 
real-time renormalization group method in frequency space (RTRG-FS^) to allow the calculation 
of dynamical correlation functions of arbitrary dot operators in systems describing spin and/or or- 
bital fluctuations. The resulting two-loop RG equations are analytically solved in the weak-coupling 
regime. This implies that the method can be applied provided either the voltage V through the dot 
or the external magnetic field /iq are sufficiently large, max{l/, ho} 3> Tk, where the Kondo temper- 
ature Tk is the scale where the system enters the strong-coupling regime. Explicitly, we calculate 
the longitudinal and transverse spin-spin correlation and response functions as well as the resulting 
fluctuation-dissipation ratios. The correlation functions in real-frequency space can be calculated 
in Matsubara space without the need of any analytical continuation. We obtain analytic results 
for the line-shape, the small- and large-frequency limits and several other features like the height 
and width of the peak in the transverse susceptibility at f2 « ^, where h denotes the renormalized 
magnetic field. Furthermore, we discuss how the developed method can be generalized to calculate 
dynamical correlation functions of other operators involving reservoir degrees of freedom as well. 

PACS numbers: 05.10.Cc, 73.63.Kv, 75.40.Gb 



I. INTRODUCTION 



The single-impurity Kondo modelS is unquestionably 
one of the most important models studied in condensed 
matter physics over the past decades. The investiga- 
tion of its equilibrium properties has caused the devel- 
opment of important theoretical tools such as renor- 
malization group methods or the Bethe Ansatz for im- 
purity systems^iii^i^. More than 20 years ago it was 
also realized^i^'^ that the Kondo model can be used to 
describe transport experiments through quantum dots. 
The developments in the ability to engineer devices on 
the nanoscale has led to the experimental realization 
of Kondo physics in such system o"'^°i^^'^^i^^'^^i^^'^^ . One 
particular advantage of these quantum dots is the al- 
most full control over system parameters like temper- 
ature, bias and gate voltages, magnetic field, and ex- 
change couplings. These possibilities have triggered 
a great interest in the theoretical study of quantum 
dots out of equilibrium. A wide range of theoretical 
methods has been applied in the past, including non- 
equilibrium perturbation theorjiiii^iiSi2ai^i2^iS^i2ii2^i2£, 
the flow-equation method22i2^i22ii2£. Coulomb gas rep- 
resentations^ii^, real-time and functional renormaliza- 
tion group methods^ '^'^i'^^'^^i'^^'^^'^^ , a non-equilibrium 
extension of the numerical renormalization group (NRG) 
method^i^, and time-dependent density matrix renor- 
malization group (DMRG) techniques^i^i^i^'^. These 
studies established the importance of relaxation and 
decoherence effects for the understanding of non- 
equilibrium physics. From an experimental point of view 
various measurable quantities like the steady-state cur- 
rent, the magnetization, and the static susceptibility have 
been calculated. 



Beside the application of these techniques there were 
also attempts to employ the known integrability of cer- 
tain impurity models in equilibrium, notably the Ander- 
son impurity model and the interacting resonant level 
model, to investigate their non-equilibrium properties. 
Konik et al.— »^ calculated the differential conductance 
in the Anderson impurity model by combining the well- 
known scattering states of the equilibrium system^ with 
a Landauer-Biittiker formalism. In this work the chem- 
ical potentials in the leads were coupled to dressed ex- 
citations rather than free electrons and the calculation 
was restricted to a subset of the excitations. In con- 
trast, a different approach was recently put forward by 
Metha and Andrei^ to treat the interacting resonant 
level model. They constructed a new set of scattering 
states of Bethe-Ansatz form which share the quantum 
numbers of free electrons in the incoming channel, hence 
allowing the application of a finite voltage in the usual 
manner. However, questions concerning the existence of 
these scattering states and issues related to the used reg- 
ularization scheme of the theory remain open. Neverthe- 
less, the interacting resonant level model has become one 
of the benchmark systems in the study of non-equilibrium 

p^..,^^j, „42.49.50.51.52.53.54.55.56 

Despite the large number of studies of impurity mod- 
els out of equilibrium only few results are known for 
the dynamical correlation functions. The spin dynam- 
ics of a non-equilibrium quantum dot has been studied 
by using a Majorana fermion representation'^"'^^, which 
yields the qualitative low-frequency properties of the cor- 
relation functions. The transverse susceptibility in a 
Kondo model was studied by Paaske et ali^ using non- 
equilibrium perturbation theory together with a pseudo- 
fermion representation of the Kondo spin. They showed 
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that the Fourier transform of the transverse susceptibil- 
ity possesses a peak if its frequency equals the value 
of the applied magnetic field, ~ ho, and that the 
width of this peak is given by the transverse spin re- 
laxation rate Their derivation was, however, re- 
stricted to either the regime Hq ^ r2 or max{|ri — 
/io|,r2} ^ ma.x{ho,V}, where V denotes the apphed 
voltage. Very recently, Fritsch and Kehrein^^''^^ applied 
the flow-equation method to study the longitudinal cor- 
relation function as well as the magnetization and T- 
matrix in a Kondo model in and out of equilibrium. The 
numerical solution of the two-loop scaling equations al- 
lowed them to study the correlation function for all com- 
binations of the parameters voltage, temperature, and 
magnetic field, provided the weak-coupling condition (i.e. 
the presence of a large enough infrared cutoff) was sat- 
isfied. In particular, these numerical solutions were used 
to compare the effects of an applied voltage and a finite 
temperature, revealing qualitative differences such as the 
appearance of Kondo splitting in the non-equilibrium sit- 
uation. In general, however, this method cannot provide 
analytic expressions for the line shape. 

In this article we will generalize the real-time renor- 
malization group method in frequency spaced to allow 
the calculation of dynamical correlation functions of arbi- 
trary dot operators in systems describing spin and/or or- 
bital fluctuations. In this setting the quantum dot is cou- 
pled to non-interacting leads which are held at different 
chemical potentials. The derived two-loop RG equations 
can be solved analytically in the weak-coupling regime. 
Explicitly, we calculate the longitudinal and transverse 
spin-spin correlation and response functions in a two- 
lead Kondo model in a magnetic fleld /iq up to order J^. 
Here Jc denotes the effective coupling at the energy scale 
Ac ~ max{y, ho} where the flow of the couphng constant 
is cut-off. In order to satisfy the weak-coupling condition 
Jc ^ 1 either the applied voltage or the magnetic field 
have to be sufficiently large compared to the Kondo tem- 
perature, Ac ^ Tk, where Tk is the scale where the 
system enters the strong-coupling regime. We note that 
the applied formalism does not rely on a fermionic rep- 
resentation of the Kondo spin but rather deals with its 
matrix representation in Liouville space directly. The 
longitudinal response function possesses a peak at the 
spin relaxation rate Fi, which gets suppressed in a finite 
magnetic field. Interestingly, in the case of a strong mag- 
netic field, V < h where h = {1 — Jc + . . ■)ho denotes 
the renormalized magnetic field (see Eq. p74p for the 
precise value), the longitudinal correlation and response 
function show "kink-like" structures at the frequencies 
fl = h^h ±V, which were also observed using the flow- 
equation method^. Here we additionally provide the line 
shape close to these "kinks" and show that the real part 
of the response functions shows characteristic logarithmic 
features a,tU = h, hztV. Furthermore we study the longi- 
tudinal and transverse fluctuation-dissipation ratios. As 
expected these ratios show a revival of the fluctuation- 
dissipation theore m^^i^^ provided the applied voltage is 



small compared to the magnetic field or the considered 
oscillation frequency. 

This article is organized as follows. In the next two 
sections we will define the general set up we want to 
study and define the used notations. This will include 
the Kondo model, the notion of Liouville operators as 
well as the definition of the symmetrized correlation func- 
tion and susceptibility. In Sec. IIVI we will then derive 
perturbative expansions for the kernels needed to calcu- 
late these correlation functions. This will be performed 
in Liouville space; the expansion is done in powers of 
the exchange coupling between the dot and the reser- 
voirs. These perturbative expansions can be applied to 
any model describing spin and/or orbital fiuctuations as 
well as correlation functions of arbitrary operators. In 
Sec. |V] we will use these results to derive the RG equa- 
tions for the kernels of pure dot operators. In the fol- 
lowing section we will further specialize to the two-lead 
Kondo model, where we use the explicit expressions of the 
Liouville operator and the coupling between the dot and 
the reservoirs to derive analytic results for the effective 
kernels appearing in the correlation functions. Finally, 
these expressions for the kernels are used in Sees. IVIII 
and lVIIll to calculate the longitudinal and transverse cor- 
relation and response functions. 

II. KONDO MODEL 

The real-time renormalization group in frequency 
space was applied in Ref. [s^ to calculate various quan- 
tities including the spin relaxation and dephasing rates, 
the renormalized magnetic field, the magnetization and 
the current in the anisotropic Kondo model in a finite 
magnetic field out of equilibrium. In this reference all 
notations which we will use in the following were origi- 
nally set up. In order to increase the readability of the 
present manuscript we will briefly recall the basic formu- 
las and notations. 

We consider a quantum dot with fixed charge in the 
Coulomb blockade regime coupled to external reservoirs. 
As shown in detail in Ref. 37, a standard Schrieffer- Wolff 
transformation leads to a Hamiltonian of the form 

H = Hres + Hs + V = Ho + V, (1) 

where Hres is the reservoir part, Hs characterizes the iso- 
lated quantum dot, and V describes the coupling between 
reservoirs and quantum dot. They are given explicitly by 

Hres = ^ J did (u) + fj,a)a+^{u})a^^{uj), (2) 

}y~a.a... 

Hs - (3) 

s 

V = J duj J dui' g,ji,^^,t,,{Lj,Lj') 

rjr}' vv' 

X : a,,^(w)a,,v(w') : . (4) 
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Here, the fermionic creation (?? = +) and anni- 

hilation (ji = —) operators in the reservoirs and v is an 
index characterizing all quantum numbers of the reservoir 
states, which contains the reservoir index a = L, R = ± 
and the spin quantum number a =|,|= ±. We mea- 
sure the energy uj of the reservoir states relative to the 
chemical potential fia of reservoir a. The eigenstates 
and eigenenergies of the isolated quantum dot are de- 
noted by \s) and Eg. The interaction V is quadratic in 
the reservoir field operators, which arises from second 
order processes of one electron hopping on and off the 
quantum dot coherently. This keeps the charge fixed and 
allows only spin/orbital fluctuations. The coupling ver- 
tex g^i,^,,'^' (w, w') is an arbitrary operator acting on the 
dot states. It is written in its most general form, depend- 
ing on the quantum numbers and energies of the reservoir 
states in an arbitrary way. As explained in Ref. [Ij, the 
RG approach can be set up in its most convenient form 
if one assumes that the frequency dependence of the ini- 
tial vertices is rather weak and varies on the scale of the 
band width D of the reservoirs. For the model we have 
in mind, the isotropic spin-i Kondo model in a magnetic 
field, this is certainly satisfied. Therefore, we will assume 
this in the following and introduce below [see Eq. ([T^ ] a 
convenient cutoff function into the free reservoir Green's 
functions. 

To achieve a more compact notation for all indices, we 
write 1 = rjiyuj and sum (integrate) implicitly over all 
indices (frequencies). The interaction is then written in 
the compact form 



V = igii' : aiay 



(5) 



: • • • : denotes normal-ordering of the reservoir field oper- 
ators, meaning that no contraction is allowed between 
reservoir field operators within the normal-ordering. 
Within the normal-ordering of Eq. ^ , the field operators 
can be arranged in an arbitrary way (up to a fermionic 
sign) , therefore the coupling vertex can always be chosen 
such that antisymmetry holds: 



511' = -ffi'i- 



(6) 



Furthermore, due to the hermiticity of V, the vertex has 
the property 



5n' = .9i'i, 



(7) 



where 1 = — ry, i/, u. 

The specific model we want to study is the isotropic 
Kondo model in an external magnetic field ho > (see 
Fig. H]). In this case the above relations read explicitly 



Hs = h„S', (8) 

1 J (Jq«')oS'V^^, for ri = -ri' = 

2 1 -(Ja'a)o^V;,^ for rj = -rj' = 



911 



, (9) 



where i S {x,y,z}, 5'* is the i-component of the spin-i 
operator of the quantum dot, cr* is a Pauli matrix, and 
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S=l/2 






FIG. 1: Isotropic spin-1/2 Kondo model coupled via exchange 
to two reservoirs. Jll and Jru involve exchange between 
the electron spins of the left/right reservoir and the local 
spin, Jul = Jlr transfers an electron from one reservoir to 
the other during the exchange process. We assume that the 
Kondo model was derived from an Anderson impurity model 
via a Schrieffer- Wolff transformation, which implies the rela- 
tion Jll Jrr ~ Jrl- 



{Jaa')o are the initial exchange couplings. We will be 
interested in the antiferromagnetic model here, i.e., we 
assume {Jaa')o > initially. If one derives the Kondo 
model via a Schrieffer- Wolff transformation from an An- 
derson impurity model (see, e.g., Ref. i37il . one further 
finds 



■Jo, 



E 



(10) 



Although the general formalism and many of the follow- 
ing formulas are also valid for an arbitrary number of 
reservoirs, we will consider the case of two reservoirs only 
with chemical potentials given by 



V 



V 
2"' 



(11) 



where V is the applied voltage which we assume to be 
positive, V > 0. 

A contraction is defined with respect to a grand- 
canonical distribution of the reservoirs, given by 

aiav = (aiai')p„a = (Jn/ /9(w) /„(?7cj). (12) 

fa{uj) = (e'^/^° + 1)^1 = 1 - fa{-oj) is the Fermi 
distribution function corresponding to temperature Tq, 
(note that the chemical potential does not enter this for- 
mula since lo is measured relative to fia)- Furthermore, 
6ii> = 5nri'5^u'5{uj — lu') is the (5-function in compact no- 
tation. Furthermore, we have introduced the cutoff by 
the band width D into the reservoir contraction via the 
density of states 



D2 



(13) 



In order to calculate the dynamical spin-spin correla- 
tion functions we have to know the time evolution of the 
density matrix p(t) . Formally, this follows from the solu- 
tion of the von Neumann equation 



Pit) 



-iH{t-to) 



JH{t-to) 
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-iL(t-to) 



(14) 



III. CORRELATION FUNCTIONS 



where 



[H,. 



(15) 



is the LiouviUian acting on usual operators in Hilbcrt 
space via the commutator. The form ([1]) of the Hamilto- 
nian yields a similar decomposition of the LiouviUian, 



L = Lr 



(16) 



with Lres = [Hres,-]^, L^s = [Hs,-]-, and Ly 
[Hv, ■]_■ We would like to note that the concept of Li- 
ouville space and superoperators have been used in vari- 
ous contexts, for example in quantum statistical mechan- 



ic 



.59.60 



Initially, we assume that the density matrix is a prod- 
uct of an arbitrary dot part ps{to) and a grandcanonical 
distribution p^es for the reservoirs, 

Pito) = Ps{to) Pres- (17) 

Furthermore, we introduce the Laplace transform 
dte''-^*-'°^ p{t) = - ' 



■ Pito) 



(18) 



where we will frequently use the notation z 
The stationary density matrix is defined as 

p"* = lim p{t) = lim p{t), 



E 



(19) 



which is understood in the sense Tr(Op''*) — 
limt^oo Tr(Op(i)) for any local operator O, and can be 
calculated using 



-i lim z p(z) — lim 

z^iO+ z^iO+ Z — L 



Pito) 



(20) 



The existence of a stationary state was proven in Ref. l2a 
using non-equilibrium pertubation theory to all orders as 
well as in Ref. [l] using the RTRG-FS, which in particular 
clarified the generation of the relaxation and dephasing 
rates under the RG flow. The reduced density matrix of 
the dot is obtained by tracing out the reservoir degrees 
of freedom 



psiz) = Trres Piz) 



-eff 



[z) 



Psito), 



(21) 



where ^^^{z) denotes the effective LiouviUian of the 
quantum dot formally defined in (|57p below. The station- 
ary reduced density matrix can then be obtained similar 

to (uni), 



Ps 



- lim ps{t) 



lim — 

z-,iO+ z - 



-eff 



iz) 



Psito)- (22) 



The quantities of interest in this article are the two- 
point correlation function of two operators A and B — 
as well as their dynamical susceptibility with respect to 
the steady state, 

SABit) = \{[A{t)^^{A)^^^B{Q)^-{B)J^)^, (23) 

XABit) = ie(t)([A(t)H,i?(0)H]„)^^. (24) 

where 



lim Tr 

t{) — > — oo 



(Oe'^*V(io) 



lim Tr Op(0) 



(25) 

Here the trace is taken over the dot states as well as the 
reservoir degrees of freedom, Tr = TrgTr^es. The time- 
evolution of the operators in the Heisenberg picture is 
given by 



A{t)n^e'"'A 



-iHt 



A. 



(26) 



Instead of calculating (|23p and ((24|) in real time we will 
study their respective Fourier transforms 



SABi^) 
XABi^) 



dte'^'^SABit), 

) 

dte'^' XABit), 



(27) 
(28) 



where — il ± iS ior t > {t < 0). The susceptibility 
admits the standard decomposition XABi^) — XabI^) + 

-^x:'ABi^). 

In order to calculate Sab^S^) and xab^S^ we introduce 
the auxiliary correlation functions 



Cl^i^)= \ die-'^"([A(0)H,S(OH]±) , (29) 



with = fi-l-i^. Its relations to the correlation functions 
are given by (see App. E| 

SABi^) - ReC\s(Sl)~2i,{A)^lB)j{^\ (30) 
XABi^) = -^C^ABi^). (31) 

The static susceptibility is related to the dynamical sus- 
ceptibility via 



(32) 



where M = (S'^)st denotes the magnetization. 

Some general properties of the correlation functions 
can be obtained by considering their spectral representa- 
tions. Let {\n)} be a complete set of basis states of the 
full Hamiltonian H, i.e. H \n) = En\n). Furthermore, 
the stationary density matrix p"* satisfies [H,p^*']_ = 0, 
thus the basis states In) can be chosen such that 



{n\p^'\m)=p:lS„ 



(33) 
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Using this one easily verifies the spectral representations 
C^Bi.^) = Y.(Pn±Pt) {n\A\m) {m\B\n) 

■ (34) 



E„ 



These relations imply 5^^t(f^) > as well as Sba{^) — 
SAB{-n), x'ba{^)^ =^ Xab(-^^), and x'ba{^) = 
~Xa_b(~^)- -'^'^ equilibrium the matrix elements of the 
density matrix are given by p„ = e^^"/^/Z with the par- 
tition sum Z, which implies the well-known fluctuation- 
dissipation thcorem^i^ 



xV(f^)=tanh^5^^t(n) 



(35) 



as well as Ux'aa'^-{^) > 0- 



IV. PERTURBATIVE EXPANSION FOR THE 
CORRELATION FUNCTIONS 

In this section we derive a perturbative expansion in 
Liouville space for the auxiliary correlation functions 
C^^(f)), which wiU serve as the starting point for the 
derivation of the RG equations below. A similar per- 
turbative expansion for the effective Liouvillian of the 
quantum dot L'g^ has been derived in Refs. [lIlH. We 
will generalize these results to Cj^(ri) while closely fol- 
lowing the presentation of Ref. 38. 

As starting point to set up the formalism we assume 
that the operators A and B admit a representation sim- 
ilar to (O, 



1 

— ai., 
to! 



: ai 



ai . . 



(36) 

where we recall the short-hand notation 1 = rjiyuj and 
sum (integrate) implicitly over all indices (frequencies). 
We further assume the operators A and B to be bosonic 
which implies m and n to be even. Eventually we will 
be concerned with the correlation functions of the spin 
operators on the dot, i.e. A,B = S~^, S~ , S^. In this case 
the operators do not couple dot and reservoir degrees of 
freedom and hence only the terms with to = n = are 
non- vanishing. However, we will keep the general forms 
pop throughout this section, which for example include 
the case of current operators^''^^ where m = n = 2. 

In order to set up the perturbative expansions in Li- 
ouville space we define the operators 



La^-[A,.]^, L^ = i[B,.]^ 



(37) 



Then using (|26p together with p4)) we obtain after some 
algebra 



Cji3(^)=(-i) lim / dte- 



X Tr 



(38) 



In the next step we use ([T9| and (|20|) and furthermore 
perform the Laplace transform t ^ D, in (j38p to obtain 



-i hm Tr La „ 



■Pito] 



(39) 

Here we have used = $7 -I- i^ to ensure convergence of 
the integral. The limit ^ iO+ has to be taken before 
(5 — !■ in order to reach the stationary state. 

The next step is to expand the expression (|39p in the 
interacting part Ly of the Liouvillian and to integrate 
out the reservoir part. This procedure was outlined for 
the reduced density matrix ps{z) in detail in Ref. Jj; we 
will generalize this to the case of the auxiliary correlation 
functions ([39|) here. First, using L = Lq + Ly with Lq = 

Era 



^5 ^ , ()39p can be formally expanded in Ly , 



C±^(f7) = -i hm^^Tr 



La t: ;— ( Ly 



r2 — Lq \ — Lq 



X Li 



Li 



1 



(40) 



Second, in order to integrate out the reservoir degrees of 
freedom we write Ly in the form 



(41) 



where we implicitly sum (integrate) over 1 = rji/uj as 
well as = ±. is a quantum field superoperator 
in Liouville space for the reservoirs, defined by (C is an 
arbitrary reservoir operator) 



ai C for p = + 
C fli for p = — 



(42) 



Here p = ± serves as an auxiliary index which is similar 
to the Keldysh index indicating whether the field opera- 
tor is acting on the upper or the lower part of the Keldysh 
contour. G^y is a superoperator acting in Liouville space 
of the quantum dot, and is defined by (C is an arbitrary 
operator on the quantum dot) 



r'PP' r - A / 911' G for p 
u^y o - Opp, I _^cgiy forp 



(43) 



In the same way we define 
1 



L, 



L 



to! 
1 



,^P....P„,_4P..^P™ : Jf-.-J^":, (44) 
^P....P„ (g^)Pi..^. : JfV..JP'.:, (45) 



where the dot superoperators Ay ''^"" and {B±)^^"^'^ act 
on arbitrary dot operators G as 



ay.. 
G ay .. 



PlPm 

n G, Pi = + 

Pi = - 



(46) 
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FIG. 2: (color online) Example of a diagram contributing 
to the auxiliary correlation function Cj[g{0,). The time di- 
rection is to the left. Each vertex G is represented by two 
adjacent black dots indicating the two reservoir field opera- 
tors associated with each vertex. The vertices A and B are 
represented in the same way by red dots. For this example 
we have chosen m = n = 2 in (|36|l as it would be the case for 
a current-current correlation function. The black horizontal 
lines connecting the vertices denote the free time propagation 
of the quantum system, leading to the resolvents ^ ^^^^ 

in Laplace space. The green lines are the reservoir contrac- 
tions arising from the application of Wick's theorem. The 
vertical blue lines between the vertices are auxiliary lines to 
determine the energy argument Xi of the resolvents. 



i<5. 



PlP2 ■ ■ ■ "PlPn 

bl..,nC, Pi 



X 



±Chi,,,n, Pi = 



(47) 



For TO = or n = we define AC — ^ [a, C]^ and 
B±C = C]^, respectively. As we consider only 
bosonic operators A and B which change the number 
of fermions by an even integer, the sign-superoperator is 
given by 



P2Pi---Pn 



(48) 



This operator was introduced to compensate additional 
signs due to interchanges of fcrmionic reservoir field op- 
erators as explained in detail in Ref. 1. 

Inserting the representations as well as ([33]) and 
(|^5)) into ([301) and shifting all reservoir field superopera- 
tors Jf to the right using 



(49) 



where we have introduced the short-hand notation Xi — 
rji [oji -f ^cii ) , one can showi that each term of perturba- 
tion theory can be written as a product of a dot part and 
an average over a sequence of field superoperators of the 
reservoirs with respect to pres ■ Evaluating the latter with 
the help of Wick's theorem, one can represent each term 
of the Wick decomposition by a diagram (see Fig. [2] for 
an example) describing a certain process contributing to 
the auxiliary correlation function (17) . Each process 

consists of a sequence of interaction vertices G^, between 
the dot and the reservoirs, and a free time propagation 
of the dot in between (leading to resolvents in Laplace 
space). Since the reservoirs have been integrated out, 
the vertices are connected by reservoir contractions (the 
green lines in Fig. [5]). This means that the various dia- 
grams represent terms for the effective time evolution of 
the dot in the presence of dissipative reservoirs. Each di- 



agram for the auxiliary correlation function has the form 



A- 



n 



7 I lim £ 



1 



. G- 



-(0) 



1 



1 



where L 



(0) 



(0) 



G- 



X2 

1 



(0) 



G 



B± 



(0) 



'S 



(50) 



G = G 



'PiPi 



indicates an interac- 



reservoir field superoperators, defined by 



tion vertex, and 7 = 7^'''^ is a contraction between the 



'11' 



— 1 1' — y ^^res 'J \ 'J \i pre 
= SlV P{i^)p' faiVP'^)- 



(51) 



We stress that only the initial reduced density matrix of 
the dot ps{to) defined in ([T7)) appears in ([50)1 as we have 
already performed the trace over the reservoir degrees of 
freedom (and hence Pres) to obtain the reservoir contrac- 
tions 7. To factorize the Wick decomposition, a fermionic 
sign has to be assigned to each permutation of reservoir 
field superoperators, indicated by the sign factor (—1)^" 
in ([50]) . For each pair of vertices connected by two reser- 
voir lines, a combinatorial factor i occurs, leading to the 
prefactor ^ in ([50)) . The value of the frequencies Xi in 
the resolvents between the interaction vertices is deter- 
mined by the sum over all variables x = r]{uj+fia) of those 
indices belonging to the reservoir lines which are crossed 
by a vertical line at the position of the resolvent (the blue 
lines in Fig. Thereby, the index of the left vertex has 
to be taken of the corresponding reservoir line. For exam- 
ple, the diagram shown in Fig. [51 is given by (the obvious 
dependence on the Keldysh indices has been omitted for 



simplicity, i.e. 7^ 



7f^*^' and Gij 



•PiPj 



r~iP-i 



i lim C(7i6723745 A2ni2(r2) 6*341114(11) G56 

f— >iO-l- \ 



fl - Lf V 2 



7:77,10789 ^±,78 Hts (£) G 



9,10 



PS {to), 



where the resolvents are defined by 

1 



ni...„(z) 



Zl...n + ^l...n — L 



(0)-' 



with 



as well as 



Zl...n = -2 + /ij 
i=l 



wi...„ = ^ Wi, Pt=riiPai, u)i = riiUi 

i=l 



(52) 



(53) 



(54) 



(55) 




7 8 9 10 11 12 

S±,78 



FIG. 3; (color online) Example of a diagram contributing to 

ED. 




FIG. 4: (color online) Example of a diagram contributing to 



As can be seen from the example (|52|) , each diagram con- 
sists of a sequence of irreducible blocks (where a vertical 
line always cuts at least one reservoir line) and free re- 
solvents l/i^l — ^5'') or l/(^ — L5'') in between. Now 
there arc two possibilities: (i) The vertices A and B do 
not belong to the same block (see Fig. [3] for an exam- 
ple), (ii) The vertices A and B belong to the same block 
(see Fig. 0] for an example). In the first case (i) one 
can formally resum those terms between the vertices A 
and B which are not connected to them similar to Dyson 
equations with the result 



(56) 



where 



-eff 



(0) 



S(f7). 



(57) 



Here the kernel I](il) contains the sum over all irreducible 
diagrams, 



m) ^ i(-i)"'(n7j.„ 

1 



xG 



G ...G- 



1 



Q + Xr- L 



(0) 



G, (58) 



where the subindex irr indicates that only irreducible 
diagrams are allowed where any vertical line between the 
vertices cuts through at least one reservoir contraction. 
We further introduce irreducible blocks T,a{^) as well as 
S^(r2,^) which are given as the sum over all irreducible 
diagrams containing the vertices A and S±, 



n 



1 



xG- 



(-1)^^ (n^ 



1 



n + Xi 
1 



L 



(0) 



G ... G- 



n + Xr- L 



(0) 



B± 



G ...G- 



1 



X,-L 



(5)G. (60) 



Obviously, in the case of spin operators, A,B G 
5" , S"^}, the vertex A has no external legs and hence 
S^(ri) = A. In contrast, although the vertex B± does 
not possess external legs either, there exist irreducible 
diagrams containing B± and at least one vertex G to the 
left and one to the right of B±. If we now proceed by 
resumming the irreducible blocks right to (and not con- 
nected to) the vertex B± similar to ([56|) and perform the 
limit ^ \Q+ using \22\. we deduce that all terms of 
type (i) contribute to 



-iTr< 



1 



5]|(f7,iO+)p|* 



(61) 



In the second case (ii) we introduce a kernel similar 
to 159]) and ((60)) which contains all irreducible diagrams 
containing both vertices A and S, 

■ - ^(-1)^'' (n- 



7 



xA- 



1 



1 



(0) 



G ... G- 



1 



e + Xr+l - L 



(0) 



G ... G- 



n + Xr- L 
1 



(0) 



B± 



C + Xs-L 



(0) 



G. (62) 



In the case of spin operators, A,B e S*^, S*^}, there 
exist no irreducible diagrams connecting A and B± , hence 
E^^(17, ^) = in this case. Now using again (|22|) for the 
sum of the irreducible blocks right to (and not connected 
to) the vertex B± we deduce that all terms of type (ii) 
contribute to 



i Tr.c 



MS 



(63) 



Hence, taking together (i) and (ii) we finally arrive at the 
main result of this section, 



-iTrs 



S|(f],iO+)p|* 



s±^(r!,io+)p|* 



(64) 



where the kernels are defined by ([58|) . ([59)) . ((60)) . and 
([S^ . respectively. 

In addition we note that the diagrammatic series can 
be partially resummed by taking all closed sub-diagrams 
between two fixed vertices together which contain only 
contractions connecting vertices between the two fixed 
ones. This has the effect that the resolvents in (|59l) . 
(|60)) . and ([62)) . are replaced by 



1 



1 



8 



(and similar for il i.e., the full effective Liouville 

operator occurs in the denominator. In this formulation 
the number of diagrams is reduced, i.e. diagrams con- 
taining closed sub-diagrams between two vertices are no 
longer allowed. 



V. GENERIC RG EQUATIONS 



In this section we will set up the generic RG equations 
for the kernels I]^(ri), I]^(il, ^), and Si^^(il, ^), for spin 
operators on the dot in a model with spin/orbital fluc- 
tuations. Hence we will assume the form (|41l) for the 
coupling between the reservoirs and the quantum dot, 
but will keep the vertex Gff, arbitrary at this stage. The 
derivation will require some relations between the initial 
Liouvillian L^s \ the vertex B± and the effective dot Li- 
ouvillian L'^''^ which can be explicitly checked for the 
Kondo model to be studied in the next section but have 
to be assumed here. These relations are (|112l) . (I132p . 
P33p . and (|B4|) . The generic RG equations for the ver- 
tex G and the effective Liouvillian L'g^ have been derived 
and solved in Ref. [ssl . we will quote these results without 
derivation when they are needed. 

Furthermore, we will restrict ourselves to the calcula- 
tion of the dynamical spin-spin correlations only, i.e. we 
will assume A,B^ {S^ , , S'^} in what follows. This 
implies in particular, that the initial values of the ver- 
tices A and B± defined in and do not possess 
any external lines, i.e. m = n = 0. Explicitly, 



-4 = - [A, 
B± = i[B,.] 



± • 



(66) 
(67) 



We will see below that this form of the vertex A is con- 
served under the RG flow. In contrast, a new effective 
i3-type vertex B±^iii with two external lines will be gen- 
erated. The fact that the initial vertex A has no external 
lines directly implies the final results for the kernels ([59]) 
and (|62p. namely 



= A, 

= 0. 



(68) 
(69) 



Hence, in the following we have to consider the kernel 
([50)1 only. We note that the results of this section re- 
main valid for any pure dot operators A and B, i.e. any 
operators ((36)) with m = n = 0. 

The RG procedure is divided into two steps. In the first 
step we will integrate out the symmetric part of the reser- 
voir contractions jff, . The reason for this is as follows: 
The effective dot Liouvillian L'^^ {z) can be diagonalized 



L%^^{z)^Y.^,{z)P,{z), 



(70) 



where \i{z) and Pi{z) denote the eigenvalues and corre- 
sponding projectors, respectively. This diagonalization 
implies for the resolvents 



1 



L%ff{z) 



E 



1 



z - \i{z) 



P^{Z). 



(71) 



Now there exists a zero eigenvalue Ao(z) — whose eigen- 
state for z iO-l- corresponds to the stationary state. 
The appearance of this zero eigenvalue can lead to in- 
frared divergencies of the frequency integrations in the 
perturbative expansions for the vertex G and the effective 
Liouvillian Vg^ as is elaborated on in detail in Ref. [l|. 
However, after the discrete RG step we can trivially sum 
over the Keldysh indices p and p' by introducing 



Gil' — E^ii'' ^11' ~ P ^11' ■ 



(72) 



In the resulting RG equations only the symmetric vertex 
G will appear. This vertex has the important property 



Po(^)Gii' =0, 



(73) 



which is independent of the model specifics. Hence, af- 
ter the discrete RG step the zero eigenvalue can no longer 
appear in any resolvent standing left to G (i.e. in no resol- 
vent except the one standing left to the vertex B±). This 
resolves the problem of infrared divergent internal fre- 
quency integrations, as the remaining eigenvalues Xi{z) 
have a strictly negative imaginary part (see l|167p and 
(|168p ). We would like to refer to Ref. [l| for a general 
discussion of this topic. 

In the second step we introduce a cut-off A into the 
reservoir contractions via the Fermi function. We then 
integrate out the reservoirs by sending A — > 0, which re- 
sults in a description of the system in terms of effective 
dot quantities like L'g^ . This second, continuous, RG 
step is further divided into two substeps, first we inte- 
grate out the reservoir degrees of freedom down to an 
energy scale Ac, and second we complete the flow from 
Ap down to A = 0. 



A. Discrete RG step 

In the first discrete RG step we integrate out the sym- 
metric part |(/a(w)-l-/Q-(— w)) = \ of the Fermi function 
in the contraction ([5T|) . The discrete RG step for the ker- 
nel ([55)1 and the vertex G has been performed in Ref. [H. 
Here we will derive the analog results for the kernel ([60)) 
and the vertex B±^iii. This is achieved by decomposing 
the contraction (jST)) according to 



pp 
liv 



5ivp ll + ^ii'7i> 



7i' = i^Pi^), li = P{^) 



(74) 
, (75) 



with uj = rjuj. Using this decomposition in ()60[) . one 
finds that each diagram decomposes into a series of blocks 
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VB± I'l 



B± 2 1' 



FIG. 5: (color online) The lowest order diagrams for the kernel 



Eg'"(r2, ^) (left and middle diagram) and the effective vertex 
3%. 1]^/ (fi, 1^) (right diagram) when the symmetric part of the 
contraction is integrated out. a (a) denotes the symmetric 
(antisymmetric) contraction 7* (7")- 



which are irreducible with respect to the symmetric part 
7* (i.e., any vertical line hits at least one symmetric con- 
traction) and connected to each other by antisymmetric 
contractions 7°. The blocks which are irreducible with 
respect to 7* can be formally resummed into an effective 
kernel £^'"(0,^) and a newly generated effective ver- 
tex i3J ^-^,(17,^). The lowest order diagrams are shown 
in Fig. \5\ Using the diagrammatic rules together with 
(|74| and the convention ((54)) . we obtain for the first two 
diagrams 



^ fill' 
1 



1 



(0) 



^11' — Lg 



B± 



(76) 



L 



(0) 



np p 



^11' + CJll' — 

and for the third one (including the interchange 1 ^ 1') 
1 



/ s nPP . 
P 12 '^12 



L 



(0) 



fll2 + U>12 JJg 

1 

X 

62 



B± 



(0) 



Gi^: - (1 - 1') 



(77) 



We use here the original perturbation series ([60]) so that 
the unperturbed Liouvillian L^g'^ occurs in the resolvents. 
Performing the frequency integrations and assuming the 
band-width to be large, we obtain 



+i-Gn.6±Gi,i + 0(G3,i), (78) 

Bl,,y{n,e} = OiG^j^), (79) 

where we have performed the sum over the Keldysh in- 
dices and used ([7^ . We stress at this point that the 
kernel and the renormalized vertex B±{fl, ^) without 
external lines are identical, 



(80) 



Nevertheless we will retain the distinction between the 
kernel and the vertex in the following, as the former ap- 
pears in the final formulas for the correlation functions 
(|64p , whereas the latter appears in the diagrammatic ex- 
pressions for the r.h.s. of the RG equations. We note that 
the frequency dependence in ([80)1 is generated during the 
flow as shown below. 



After integrating out the symmetric part of the Fermi 
function in this way, we obtain a new diagrammatic se- 
ries for the kernel analog to ((60)) . The Liouvillian and the 
vertices have to be replaced by the effective ones and the 
contractions between the effective vertices contain only 
the antisymmetric part 7°. Due to ([75]) there occur no 
diagrams including the new vertex B^. £,)■ Further- 

more, since the effective quantities have become energy 
dependent (also the effective vertex becomes energy 
dependent in higher order perturbation theory) , one has 
to replace 

1 



X, - L 



(0) 



(81) 

(with z = r2,^) in (j60p . Since the antisymmetric part 
of the contraction ([7^ does not depend on the Keldysh 
indices, only the effective vertex G° averaged over the 
Keldysh indices occurs in the new perturbative series. 



B. Continuous RG equations 

In the second continuous RG procedure we deal with 
the remaining antisymmetric part of the Fermi distribu- 
tion function, where in each infinitesimal step a small 
energy shell is integrated out. Instead of integrating out 
the energies on the real axis, it has turned out to be more 
efficient to integrate out the Matsubara poles of the Fermi 
distribution function on the imaginary axis^i^^. This is 
achieved by introducing a formal cutoff dependence into 
the antisymmetric part of the Fermi distribution by 



/a(^) 



1 



^^T.(A-KI), 



(82) 



where = (2n + l)TrTa are the Matsubara frequencies 
corresponding to the temperature of reservoir a, and 



6It(w) = 



OH 



27tT 



for \uj\ 
for lujl 



> ttT 
< ttT 



(83) 



is a theta function smeared by temperature. For A = 
00, (l82|) yields the full antisymmetric part /a(w) — ^ of 
the Fermi distribution. In each RG step, one reduces 
the cutoff A by dA, and integrates out the infinitesimal 

part — fa^"^^ — of the Fermi distribution. For 

example, the new effective Liouvillian at scale A — rfA 



- A-dA 



(z) = LUz)-dLUz) 



(84) 



and similarly the new effective vertices G^y'^^{z) and 

S±7iif(r2,0 as well as the kernel s|^^"''^(17, can be 
calculated technically in the same way as for the first dis- 
crete RG step. The only difference is that an infinitesimal 
small part is integrated out, so that the RG diagrams con- 

tain only one contraction involving dA-^. Furthermore, 
since the diagrams have to be irreducible with respect to 
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1 2 



2 r 



FIG. 6: (color online) RG diagram for the renormalization 
of the vertex j^j/(Sl,z) in O(G^). The slash indicates the 
contraction where the Fermi function has to be replaced by 



1 2 g 2 1 12 2 1 12 




3 1 12 




2 3 3 1 



FIG. 7: (color online) RG diagrams for the renormalization 
of the kernel E^'*(!^, z) up to O(G^). 



this part, this contraction must connect the first with the 
last vertex of the diagram. Using this procedure the RG 
equations for the dot Liouvillian Lgjz) and the vertex 
G^^, (z) have been derived in Ref. [Sa. Here we will use 
this technique to obtain the RG equations for ly (O, ^) 

as well as I]^'^(il,^). 

The diagrams contributing to the RG equations for 

B±,iv{n,o = Bi,,,{n,o and s±(r!,0 = s±'^(r!,c) 

are shown in Fig. |6] and Fig. [71 respectively. Using the 
definition 



7f = piCo)ftiCo), 
together with the convention 

1 



ni...„(z) 



Zl...n + t^l...n — Ls{zi,,,n + UJl...n) 

we obtain the following RG equations: 



(85) 



, (86) 



dA 



-^f Gi2(17)ni2(r!)S±(17i2+Wl2,Cl2+<:Sl2)ni2(C)G2i,(a2+C^12) - (1 ^ l') 



(87) 



and 



-^l2 Gum ni2(f^) B±{n,2 + ^12 + L^U) Tl^iO G^li^U + ^ll) 



dA 
dA 



72 



B±,i2m ni2(0 ^21(62 + cd,2) + Gi2m ni2(r!) ^±,21(^^12 + ^^12, 



7^^ Gi2(r!) Ui2in) [B±ini2 + oj^, 62 + 0J12) ni2(0 633(62 + ^^12) 

+ 623(1^12 + UJ12) ni3(r!) 6±(17i3 + W13, 63 + W13) 1 ni3(0 G3i(63 + tsis). (88) 



We recall here that the kernel and the vertex without 
external lines B± equal each other, see ((80|). which yields 
a closed set of RG equations. We will further show in 
App. |B] that the two-loop diagrams for i3±,ii' as well 
as the one- loop diagrams containing B±^iii itself do not 
contribute at second order in the coupling constant and 
hence can be neglected on the r.h.s. of ([57|) . 

The initial conditions of the RG equations are given 
by ^ and (HH). Since 7^^=° = 0, the solution at A = 
provides the result for the kernel 

l]|(f},0 = S±(r!,0|A=o, (89) 

from which the correlation functions can be calculated 
via ((64)1 . 

As the resolvents and the vertices on the r.h.s. of the 
RG equations are analytic functions in all frequencies Coi 
in the upper half of the complex plane, all frequency in- 
tegrations can be calculated analytically. The only poles 



contributing are the ones of the contractions and their 
derivatives given by ([55)1 with ([5^ as well as 

dA 2n \uj — iAt^ oj + iAt^ J 

Here At„ denotes the Matsubara frequency lj" which lies 
closest to the cutoff A. After performing the integration 
we find^ that, due to the presence of the cutoff function 
p{uj) — 1 r.h.s. of the RG equations gives a neg- 

ligible contribution for A ^ D. Therefore, we can start 
the RG at Ao ~ and omit the cutoff function p{oj) (the 
precise ratio between Ag and D is determined such that 
no linear terms in D in the effective Liouvillian are gen- 
eratedfii^S). As a consequence, only the Matsubara poles 
of the Fermi function in the upper half plane contribute 
and all real frequencies are simply replaced by Matsubara 
frequencies. From now on, we write the frequency depen- 
dence explicitly and define the analytic continuation of 
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the Liouvillian and the vertices in imaginary frequency 
space by 

Giv{E,uj;lji,ujv) = Gii'(S + iw)k.-*i^,,(91) 
LsiE.Lo) = Ls{E + iLu), (92) 

6±,ii.(l^ + i<5,e + iOb.-i-., (93) 
E|(f],<5,e,r) - S|(r! + i5,e + i^'), (94) 

where we keep the real and imaginary parts of the 
Laplace variable z = E + icu and the external frequen- 



cies 57 = 57 + i(5 and ^ = ^ + i^' separated from now on. 
Furthermore, w = w", = uj!^' correspond to Matsub- 
ara frequencies and the compact indices 1 and 2 on the 
l.h.s. do no longer contain the frequencies oji. With the 
definition 

the RG equations (|87|) and (|88|) in Matsubara space can 
be written as 



dA 



B±sv (^^, 5, e'; cji,ojv) =i Gi2{fl, 5; wi, At„, ) n(f7i2, (5 + uji + At„, ) S±(57i2, <5 + uji + At„, , 62, + ^1 + At„, ) 



n(a2,^' + c^i + AT„JG2i.(62,e' + ^i + AT„,;-AT„,,wi') - (1^1'), 



(96) 



and 



— S]|(r!, ,5, e') = Gi2(r!, S; At„, , ^2) n(r!i2, <5 + At„,+ C02) 6±(r!i2, <5 + At„,+ ^2, 62, + At„,+ ^^2) 



X n(Ci2, C' + At„,+ W2) (521(62, + At„,+ W2; -c^2, -At„^ ) 
+e±,i2(17, (5, C'; At„, , W2) n(a2, + At<,,+ a;2) G2i(a2, + At<,,+ ^^2; -t^2, -At„^ ) 
+Gi2(f7, 6; At„^ , ^^2) n(57i2, (5 + At„^+ ^^2) ^±^21(^^12, <5 + At„^+ t^2, C'; -^^2, -At„ J 
-i Gi2(0, 5; At„^ , UJ2) n(17i2, (5 + At„^+ UJ2) 

S±(f7l2, (5 + At„^+ CJ2, 62, e' + At„,+ t^2) n(a2, + At„^+ C^2) G23(Cl2, C' + At„,+ 

+G23(r2l2, l5 + At„^+ W2; -W2, ^^3) n(r2i3, (5 + At„^+ Wg) S±(r2i3, <5 + At„^+ Wa, Cl3, + At„,+ W3) 

n(Cl3,?' + At„^+W3)G3i(Ci3,C' + AT„^+W3;-^^3,-AT„J■ 



(97) 



r 



In these equations we implicitly sum over all indices and 
Matsubara frequencies on the r.h.s. of the RG equations 
which do not occur on the l.h.s. Only positive Matsubara 
frequencies smaller than the cutoff A are allowed and each 
sum has to be written as 



27rr„ 



n 



which reduces to an integral dio for zero temperature. 

In the next two subsections we will solve the RG equa- 
tions (|96p and ([W)) analytically in the weak coupling 
regime up to O(G^). Weak coupling is defined by the con- 
dition that the renormalized vertices Gi2{E,uj,uji,u!2) 
stay small compared to one throughout the RG flow, so 
that the expansion in powers of G on the r.h.s. of the 
RG equations is well defined. This condition is fulfilled 
if the various cutoff scales occurring in the resolvents are 
much larger than the Kondo temperature Tk at which 



the vertices would diverge in the absence of any cutoff 
scales. 



C. Weak coupling analysis above Ac 



As is discussed in detail in Refs. Illl38l there exists a 
characteristic energy scale 



Ac = max{|£;|, \fia\,h}, 



(99) 



where h » Hq is the renormalized magnetic field. For A > 
Ac the cutoff scales \E\ = \fj,a\ and h can be neglected 
in the RG equation for the vertex G (see below). This 
leads to a reference solution G^^-* which serves as the 
starting point for a systematic expansion in powers of the 
coupling constant J* , where G^j"* °c (see pB]) together 
with This yields a perturbative solution of the RG 
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equations in the regime A > Af,. These results serve 
as initial values for the flow in the second regime < 
A < Ac. Here the renormalization of the vertex G12 is at 
least of order , where = J^^^" . Provided the weak- 
coupling condition Jc ^ 1 is satisfied all quantities can 
be calculated perturbatively. This fact crucially relies on 
the appearance of some relaxation/dephasing rate in the 
resolvents ([7T|) , which is guaranteed by ([75)1 (as then in all 
resolvents standing left to a vertex G the zero eigenvalue 
of the Liouvillian cannot contribute). This analysis has 
been performed for the Liouvillian and the current kernel 
in the anisotropic Kondo model in Ref. [H. We note that 
the zero eigenvalue may appear in the resolvent left to the 
vertex B± in ([50]) . We will show below that this docs not 
lead to any problems in the calculation of the spin-spin 
correlation functions in the Kondo model up to order J|. 

Regarding the appearance of the external frequency 
as one of the cut-off parameters in (IM)l we see from the 
perturbative expansion ([60)1 that does not appear in 
the resolvents and as a cut-off parameter for all vertices 
right to B±. This fact, however, will only affect the re- 
sults in the regime ^ V,h. We will therefore use ([M]) 
as unique cut-off for all vertices appearing in the deriva- 
tion of the kernel E^. In Scc. lVII we will show that for the 
spin operator in the Kondo model the difference yields a 
correction cx 1/il and can thus be neglected. Neverthe- 
less we stress that all vertices appearing in the stationary 
reduced density matrix p|* in (|22p do not possess as 
cut-off parameter. We thus deduce that in order to stay 
in the perturbative regime we cannot rely on the external 
frequency but have to require maxjV, h} 3> Tk- 

We finally note that temperature serves as a unique 
cutoff for all terms on the r.h.s. of the RG equations as 
for A < 27rrQ the Matsubara sums are reduced to one 
term and the cutoff A^^ = ttTq becomes independent of 
A. This trivial cutoff is set to zero in the following, i.e. 
we will set = 0. 

After these preliminary remarks let us turn to the eval- 
uation of the RG equations. The one- loop RG equation 
for the vertex G is at zero temperature given bj^*^ 

■^Gii'{E,u;;u;i,ujy) = iGi2{E,Lj;LJi,A) 
aA 

X n{Ei2,uj + A + uj2)G-2riE,uj;-A,ui,) - (1 ^ 1')- 

(100) 

The RG equation for the reference solution G^^' is ob- 
tained by assuming A to be much larger than any other 
term appearing in the resolvent, which gives 



dA^"' " A 



'-'12 '-^21' ^1'2'-'21 



The initial condition for G^^^, at A = Aq ~ I? is the bare 
vertex Gw defined in The leading order solution 

is proportional to the coupling constant G*^^-* oc J (A). 
We stress that the term on the r.h.s., which is ^ J^/A, 
contributes to the change of the vertex at order J. This 
is a general feature of the RG above Ac. In order to 



(101) 



calculate the change of a quantity at order J" one has 
to analyze those terms ^ J"/A, where A ~ n,fia,h 
is some energy scale, and ~ J'^+^/A on the r.h.s. of 
the corresponding RG equation. In contrast, terms ~ 
J"+i(A/A)'=/A {k > 1) do not contribute to the change 
at order J". The RG equation for the vertex G^^'! is given 



by pOTj) with the replacement^ G^^^!,G^^^^ 



"^21" "^21 



Using this leading order solution we can formally ex- 
pand all quantities in powers of J, i.e. 



Giv{E,uj;uji,uJi>) = 

G[\]+G^^,]{E,cu;cui,Ljy) + 
LsiE,co) = 



Lf+L's^'iE,u;)+L's^>{E,u;) + 

6i',\i,(r!,5,C,e';t^i,^2) + ... 



-(1) 



(2), 



v±^(0) 



Here L^"-* — [Hs, is the bare dot Liouvillian. We recall 

dHOl), which implies i:^''-"\n,S,^,^') B';^\n,S,^,£,') in 
all orders in J. Furthermore, we have already indicated 
which terms will depend on the Matsubara frequencies, 
external frequencies + i(5 and ^ + i^' and the Laplace 
variable E + ito. 

The vertex G and the Liouvillian were calculated in 
Ref. We will here state those results needed for the 
calculation of l|104p and (jlOSp . The second order vertex 
G'^-* is further decomposed as 



(102) 
(103) 
(104) 

(105) 



G^^,]{E,uj;uji,uy) ^iGil":'^ +g[IT 



+ Gf-^^^{E,uj;uJi,u!i>). 



Here the vertex G^^"^^ is given hy^ 



•-^11' 



^12 ^21' 



'-^l'2'-^21 



(106) 



(107) 



For the Kondo model the vertex G*^^'^^-' turns out to have 
the same matrix structure as the leading order solution 
This implies that both can be put together by 
redefining G^^^ = G'^^ +G^^'^^\ which amounts to a two- 
loop renormalization of the Kondo temperature^ i^^'^^ . 
Furthermore, the vertex G^^''-' is generically given by 



G«ln: 



A 



- LOi — ii?i2 + iL 



(0) 



■G«-(l^l'). 

(108) 



A 



The zeroth-order Liouvillian is given by the initial con- 
dition Lg'' = [Hg, .]_ with ([5]) while the first-order Liou- 
villian is further decomposed as 



(109) 
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where L^^^ and Z'^^ do not depend on the Laplace vari- where we have defined 
able. The second-order Liouvillian was calculated in 
Ref. [38! , however, we will not need it for the solution 
of (HH) and (HZl) in the regime A > A^. 

Let us now turn to the calculation of (|104p and (jlOSp . 
The zeroth-ordcr term of the kernel is just given by the 
initial condition (|67p . i.e. 



/Ca(z) = In 



2A-iz 
A - iz ' 



(114) 



.±,(0) 



(0) 



(110) 



The leading term in (|113p is extracted by treating the 
terms ^ zj k separately, 



For the derivation of an RG equation for S^'''^'' we have 
to keep all terms (x J'^/A on the r.h.s. of ((97)) . Taking 
the zero-temperature limit and keeping only the order 
in the resolvents we obtain 



/Ca(z) = /Ca(2) + — , 



(115) 



dA' 



.±■(1) 



xi3. 



1 



(0) 

^ C12 + i^' + iA + icc;2 - ' 



f^i2 + i(5 + iA + icj2 - L 

Cz2i • 



(0) 



We evaluate this integral by assuming 



(111) 



where Kt^iz) can be integrated by /Ca(z) ~ ■jrFt^[z) with 



, , . , 2A-iz iz /, A(2A-iz) , , 
FA(z) = Aln— — In-)— — -/ + 1 , (116) 



K-iz 2 V 2(A - iz 



(0) r(0) _ 



= kB. 



(0) 



(112) 



which has to be checked for the specific model at hand. 
In the Kondo model we will find k = ±ho for B — 
and K = for 5 = 5'^ (see Sec. |Vll below). pT2l) can be 

used to shift S^''' to the right and evaluate the remaining 
integral by a partial fraction expansion 



V±.(i) 



j(0) 



i5-4°))Sf 



which has the asymptotic behavior Fa{z) ~ A(ln2 + 
O(zVA^)) as A ^ 00. Using (fTTS]) in (fTTS]) we obtain 
the frequency independent result 



±,(1) _ „(o) ^(1) 



(117) 



The initial condition is given by ([75)) . i.e. S^'^^''|a=Ao 
0. As mentioned above we identify this with the vertex 
in first order, B±.^ 



.±,(1) 



B^±' ICa{^i2 + ie' - Lf) G^V' (113) (HUl, which is given by 



^(1) 



Using the result for Sj.^' we derive in App.lBjthe vertex 



A + cji -i(»i2 +i^-4"^) 

12 ^ 

-(0)n 



G^V In 



G(^^)-(1<^1') 



_^(i) ^(0) 1^ A+c.i-i(a2 + ie--x, ) ^(1) + (1 ^ 



(118) 



where k = k + 0{J) (see Sec. IVII below). We note that B^2\v °^ J^/A. For large A and using the result (jllSp we can 
derive the RG equation for the kernel in second order (see App. [C| 



(0) f^' , 



-(0)^ 



(1) 



•^21 



1 

2A 



■^(2qi)„(0)^(1) ^(l)„(0)^(2ai)l 1 ^(1)^(1)^(1), 1 ^(1) 

•-^12 "± <-'2l + '-'12 "± '-'21 J^^'-'12"± '-'21 +^'-'12 



^(1)^^(0) 



G^ 



(119) 



r 



where Fj^{z) — F\{z) — i|ln2. The initial condition is decompose the second-order kernel as 
given by (|78p . At this point it turns out to be useful to 



SB''"(r!,5,e,e') = 

fl20) 
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Here T.g''-^''\n,5,^,^') is given by 



xG 



(1) 



which satisfies the initial condition 



(0). 



g: 



(1) 



^±.(2a) 



A=Ao 



-^Grl2D± G21. 



(121) 



(122) 



The remaining terms in ()120|) satisfy 



i 

~2A 

.±,(2c) 



r.(2ai) „(0) ^(1) ^(1) „(0) ^(2ai) 



(123) 



dA 



(1) 



2A 



/^(i) 



(124) 



with initial condition (according to ([75]) and (|122p ) 

^2 



.±,(26) 



A=Ao 



.±,(2c) 
^B 



A=Ao 



— Crl2 0± Lr21 



In 2 



2 G12 B^^'' G21 -I- i— G12 B'"^' G21 . 



?(0) A_ 



(125) 



The RG equations in the regime A > Ac derived in this 
section for a model describing spin/orbital fluctuations 
will be specialized and solved for the isotropic Kondo 
model in Sec. I VI Al below. In the next section we will first 
study the effect of the RG flow in the regime < A < Ac. 



In into the renormalized vertices. At A = Ac, the bare 
coupling constant is replaced by a renormalized one Jc, 
all logarithmic contributions are eliminated, and a simple 
power series in Jc remains. Thus the RG equations can 
be solved perturbatively provided Jc <C 1. In addition, 
the Liouvillian in the resolvents is replaced by the full 
effective Liouvillian L'^J^ {z). 

When calculating diagrams containing the resolvent 
()71|1 an obvious problem is that the frequency depen- 
dence of the effective Liouvillian L'g^ {z) is not known 
explicitly. To circumnavigate this complication we use 
the following approximation for the resolvents. 



1 



Here the poles Zi of the resolvent follow from the self- 
consistency equation 



Ai(2i). 



(127) 



The residue satisfy^ ai — 1 + 0{J) and hence can be set 
to one in the following. 

Starting with the one-loop RG equation ()100|) for the 
vertex G we observe that the terms on the r.h.s. are 
already of 0{J^). Therefore, the renormalization of the 
leading-order vertex G'-^-' stops at A = Ac and we have to 
use its value G^^^'^ oc Jc in all calculations from now on. 
Note that we indicate the use of the coupling constant at 
Ac by the additional superscript c. Furthermore, as we 
are eventually interested in the kernel up to 0{J^) 
we deduce that the higher-order vertices G^^''^-' and G^^**^ 
will not be needed below Ac. 

Using the replacement G^^^ G^^'>'' cx Jc in we 
obtain up to 0{J^) 



D. Weak coupling analysis below Ac 

As explained in Ref.lHthe RG above Ac has resummed 
all leading and subleading logarithmic contributions in 



E±(f7, 5, C, C) = + ^b"' + ^b" i^, S, e') 



(128) 

where we have already used that the flow of the kernel in 
order J also stops at Ac. The second-order kernel satisfies 



±,(2) 



duj2 G[\^"u{ni2,s 



A + u;2) Bf^ mi2,e + A + cc>2) G^j 



(l)c 



K.k{ ^12 + i(5 - Zj) - /CA(gi2 + iS,' - 

17-e-(z.-z,)+i(^-e) 



G^^'P.{z,)B'^^>P,{z,)G 



(0) 



(129) 



where we have used the approximation ([126]) and only obtain 
kept terms cx J^. The initial condition is given by the 

solution at A — Ac. Using the decomposition (|120[) we Tr^'^'^\Q, 5 ^ ^') 
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where the flow of E^'*'^'''' below Ac is governed by (|129p . 
The initial value at A = Ac is given by (|12ip . which we 
can rewrite as 



A=A, n-^-K + i{s-^') 

J2 [fL {^12 +iS- z,) - (62 + i?' - zj) 



In doing so we have assumed 

i 

1 = ^P,(z,) + 0(Jc), (133) 



(131) 



and neglected terms of order J^. We will show in App. |E] 
that these assumptions are fulfilled for the Kondo model 
in a magnetic field. We will further show that k = Zi — 
Zj+0{Jc) for all pairs {i, j) for which the last hne in (|13ip 
is non-zero. When applying the results of this section 
to other models of spin/orbital fluctuations one has to 
ensure the validity of the assumptions made above. The 
initial value problem (|129p with (|13ip is readily solved 
using JCa{z) = ^^a(^) with 



FA(z) = i?A(z) + | (ln^ + 1 



(134) 



The result for the kernel at A = is then obtained using 
F^^oiz) = -if In 2: 



A=0 



IE 



1 



j(0) 



(l)c 



(rii2 + i(S - z^) In 



iAc 

O12 +16 — Zi 



+ 1 -(ei2 + ie'-^,) In 



iAc 

62 + if - zj 



(135) 



Together with S^''^^-**^ and T,^'^^'^^" determined in the 
RG procedure above Ac this yields the final result for 
the kernel in second order in Jc- It is applicable to any 
operator B which does not couple the dot and reservoir 
degrees of freedom, i.e. whose initial value satisfies n = 
Furthermore the calculations were done for a 



generic model describing spin or orbital fluctuations as 
the initial vertex G was assumed to have two external 
legs. The only assumptions we have made regarding the 
model specifics are the commutation relations (I112p and 
(|B4p . Eqs. (|132p and (|133p as well as the specific relation 
between the parameter k and the poles Zi and Zj, which 
have to be determined from the self-consistency equation 

In the next section we will apply the results de- 
rived above to the spin-spin correlation functions in the 
isotropic Kondo model. In particular, we will show that 
the assumptions discussed above are justified in this 
model. Finally we note that similar results for the ef- 
fective Liouvillian have been derived in Refs.Hlii 



VI. EXPLICIT RG EQUATIONS FOR THE 
KONDO MODEL 



pling to the leads as given in ([5]) and 



A. RG flow above Ac 

The first step is to represent the initial vertex and 
Hamiltonian in Liouville space, 

rPP - I f (4a' )o LP' al^, for r; = -r/' = + 

" 2 1 -(J^Joi^V;,^ for ,7 = -V = - ' ^ ^ 

Lf = [Hs, ]_ = ho + L-^) = ho L\ (137) 

where the spin superoperators 1/ — {LP^ , LP^ , LP^) are 
defined by their action on an arbitrary operator A on the 
dot Hilbert space via 



L+A = SA, L-A 



-AS. 



(138) 



An explicit matrix representation for the superoperators 
1/ is given in App. [D] where also further superoperators 
are defined. 

The leading-order vertex G^^^ was derived in Refs.fUlBSl. 
It can be parametrized for rj = —77' = + as 



In this section we will speciahze the generic results de- 
rived above to the case of the spin-spin correlation func- 
tions in the isotropic, antiferromagnetic Kondo model 
in a magnetic field. The Hamiltonian was presented in 
Sec. ini in particular, the dot Hamiltonian and the cou- 



rt^) - ~T rr 



(139) 



where g_ is the vector formed by the Pauli matrices. The 
vertex for 77 = —77' = — is obtained using the antisymme- 
try Gil' = ^Gi'i. Inserting p39p into the RG equation 
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(jlOip and using the antisymmetry G 
(ID15p . and Jo,p = J pa we obtain 



(1) 

12 



d _ 1 



(140) 



If we assume the form pU)) . i.e. Jqq' = 2^XaXai J with 
= 1, we obtain the usual poor-man scahng equa- 
tion 



with the solution 

1 



J(A) = 



2 In^^ 



J(A)2, J(Ao) = Jo, (141) 



Tk ^ Ke-^l^^°. (142) 



Eq. (|140p explicitly shows that the term ^ J'^ / K on the 
r.h.s. contributes to the renormalization in order J. Sim- 
ilarly, the renormalization of a quantity in order J" is 
determined by the terms ~ J"/A as well as ~ J"+^/A. 
In contrast, a term J"'+^(A/A)'^/A with fc > 1 does 
not contribute at order J", as can be seen from 



, A*^ J"+i A'= '•'^ 



A"= A' k 

jn+l ^ ^ I 

~ A^ k k 
The vertex G^^! is given by 
^(1) 



dA' I — L 1 7"+i 



A A fc jn+2 



1 



(144) 



Beside the leading-order vertex (|139[) and the zero- 
order Liouvillian (|137p we also need explicit expressions 
for the vertex (5'^°^' as well as the Liouvillian in first 
order. These are given bj^i^ 



•-^11' — 'J aP 'J I3a' U 'Q-cra'^ 
1 



(1) 

= trJL'^ 



4 = ^irJhoL'^, 



(145) 

(146) 
(147) 



where J = (Jaa') is the coupling of the leading-order 
vertex G'^^' and the trace tr is taken in the reservoir in- 
dices, tr J = Jjiji + Jll- Furthermore we have taken the 
scaling limit Jq — > 0, Aq w I? ^ oo such that the Kondo 
temperature remains constant. 

As the next step it is straightforward to derive the 
following results for the initial vertex B^f''' 



(0) 

+ 

(0) 



(148) 
(149) 



where s takes the values s = z,± for S = S^,S^ and 
we have set = U^, j ~ 1,2,3. We stress that in this 



way the spin operators are directly represented by their 
matrices in Liouville space. Hence, we do not have to 
use a pseudo-fermion representation of the Kondo spin. 
The vertex A defined in ([66|) is just given hy A — h B^_^\ 



Furthermore, we recall that the kernels in zeroth order 
are just given by Syi(ri) = A and T,^{i},^) — b'"^\ re- 
spectively. Now the RG equation (I117P for S^'^^^ reads 



±,(1) 



0, 



for S 



,(1) 

B 1 



^^tr J^B^°\ forS-'W, 



(150) 



where we have used (lDT3t and ((DT6)) . The additional 
factor of two is due to the implicit summation over rj. 
Hence the solution in the scaling limit is given by 



^ Bf 



0, 



(151) 



^tr JSL"^ = -itr JL^ (152) 



In second order we will here determine only the terms 
S^'*-^''-' and E^'^^*^-* as their fiow is cut-off at A = Ac. 

The remaining term S^'^^"' will be derived in the next 
section. Hence we have to solve the RG equation (I123p . 
which using (|DT7l) and (|D18l) reads 



Ll for Et^f^^) 



dA^^ ^A"'' 1 Ll for £7'^'") 



-'B 
^B 



(153) 



Analogously we obtain for (|124p using (|147p . (|D16p and 

(IdTqI) 

:^^B''^^^=-4trJ^trj|^ f-^^':"' (154) 
dA ^ 2A I Ll for E, ^ ' 



^B 



The solutions read 



V + ,(2f') 
^B 


2 " 


^B 


- ^trj2L3 
2 


^B 


= 0, 


^B 


= iJtrJ^L^, 



(155) 
(156) 
(157) 
(158) 



where we have already taken the scaling limit so that the 
contributions from the initial condition are negligible. 



B. RG flow below A,, 

As we have already explained above the RG flow of the 
leading order solution G^^^ stops at the scale Ac which is 
given by the maximal one of the external parameters, 

Ac = max{|r2|,y,/io}. (159) 

The value of the vertex at Ac is given by (|139p with 



Jc = J(Ac) = ( "^f 



(160) 
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where due to PTT)) the couphngs satisfy 



The LiouviUian up to 0{J^) can be parametrized a; 



.38 



Jr = 2xrJc, Jl = 2xlJc, Jnd = VJrTl, (161) 



with XR + XL ~ 1 and 



Jc 



2 hi^ 



(162) 



The asymmetry ratio is defined asr = Jl/ Jr. As already 
mentioned the flow of aU vertices right to B± does not 
stop at Ac as defined in (|159p but rather at max{V, ho}. 
This affects the result for the kernel only in the regime 
n :s> V,ho. We wiU discuss the changes in this case 
separately at the end of this section. We stress, however, 
that the flow of all vertices G needed for the derivation 
of the stationary reduced density matrix p'g is cut off by 
max{V, /ig}. Hence in order to stay in the weak-coupling 
regime we need 



(163) 



As the definition of the scale Ac is to some extent arbi- 
trary as long as it remains of the order of the external 
energy scales in the problem, it is necessary to study the 
effect of a redefinition Ac — > Ac with A'c/Ac ^ 1. This 
will induce a redefinition of the coupling as 



1 



2 In 



= Jc-2j2ln^-fO(j3). (164) 



As we will show below, the redefinition Ac A^ does not 
change the final results for the LiouviUian or the kernel 
up to order J^. 



Ls{E, uj) =h{E, uj) L'^ - iT'^iE, oj) 

, , (165) 

-iT''{E,cj)L'' -ir^'{E,uj)L^'. 



The LiouviUian can be diagonalized (see (|70p ) using the 
eigenvalues {z = E + iuj) 



XoiE,Lj) = 0, (166) 

XiiE,Lo) = -ir"(£;,c^), (167) 

X±iE,Lu) = ±h{E,uj) -iT''{E,uj) -iV^iE^u), (168) 

and the projectors 



PoiE,uj) 
Pi{E,uj) 



= L"- 



T-{E,uj) 



T''{E,uj) 



(169) 
(170) 
(171) 



The eigenvalues (|166l) - (jl68p can now be used to de- 
termine the poles Zi of the resolvent defined in (|127p . 
Solving the self-consistency equation one finds zq = 0, 
zi — — i Fi — — i r° (0, 0) , and z± = ±.h — \ T2, where the 
spin relaxation and dephasing rates and the renormalized 
magnetic field are given up to 0{Jf) by 



h 



{Jl + Jl) ImH2{h) + JrJl {lmn2{V + h)+ ImH2{V - h) 



= 1- 



2 ' 2 

Jr + Jl , {Jr + Jl? 



2Imni{V) + ImT-L2{V + h) + ImT-L2{V - h) 

JrJl 



P A- T2 

K-'-^^^Ren2{h) 



(172) 
(173) 

ReH2{V + h)- Re'H2{V - h)) . (174) 



The higher order terms ~ Jc In .. . for the rates were 
obtained in Ref. 38. The voltage was defined in pT|) and 
we always assume V,ho>0. As these rates are obtained 
from the LiouviUian at z = and z = ±h, respectively, 
the external frequency Q does not appear as a cut-off in 
the definition of Ac. Furthermore, we have defined the 
auxiliary functions 



n,iE) = In 



Ac 



E^ + r? 



ii^arctan— , (175) 



which arises from terms like 



(z - Zj) In 



iAc 



(176) 



by neglecting the imaginary part of Zj, which is propor- 
tional to V-YI2 oc J^, in the prefactor and takeing the real 
and imaginary parts. We note that as z± — ±h — i r2 
the renormalized magnetic field automatically appears in 
the logarithm. We have therefore also kept h in the linear 
prefactor. The deviation of ImTii{E) from ^\E\ is only 
important for \E\ < Fj and will be neglected otherwise. 
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Furthermore, we have omitted the imaginary part of the axis. For uj = and using (|175p . the functions in the 
Laplace variable, uj, since for the correlation functions parametrization p65p are given up to order hy^ 
calculated below we only need the Liouvillian on the real 



HE) = 



1 + ^^^4^ - ^^-^^ + -^^^'^ + ^^^^ + h) - rL2{E - h) 

JrJl 



n2{E + V + h)+ H2iE -V + h)- H2iE + V-h)- 'H2{E ~V ~h) 



r'^(S) = -i(Jfl, + Jl) ( 1- 
. Jb Jl 



Jr. + Jl 



E-i 



Jl + J 



n2{E + h) + n2{E-h) 



2 J 2 
H2{E + V + h)+ n2{E -V + h)+ H2{E + V-h)+ H2{E ~V-h) 



^J^±^ I ni{E) - i-n2{E + h)- i-n2{E -h)]~i 



JrJl 



Hi{E + V)+Hi{E-V) 



~^n2iE + v + h)- ^n2{E -V + h)- ^n2iE + v -h)-^n2{E-v -h) 

^{jR + JL)^ho^^{jR + JL)^h. 



(177) 

(178) 

(179) 
(180) 



In the last line we have replaced the bare magnetic field 
by the renormalized field, such that the latter appears 
consistently in all functions (|177p - (|180p . The change is 
of O(J^). We further note that a redefinition A^, 
yields the same result for the Liouvillian with the re- 
placement Jc J'^. Naively, the linear terms in h{E) 
and T'^{E) yield additional contributions in order J'^. 
These are, however, exactly cancelled by terms appear- 
ing from the logarithms using Ht{E) = H'^iE) - Eln^, 
where 'H-(-E) is given by ifTTS)) with A^ -> A'^. The facts 
presented above allow us to justify the assumptions (|132p 
and (|133p made in Sec. IV Dl above for the specific case of 
the Kondo model (see Add. [E]). 

After the recall of the Liouvillian in second order we 
will finally evaluate the kernel in the Kondo model. 
The flow of the kernel also stops at Ac except for one 
term, namely Hence we flnd 

I]|(f7,<5,^,0-S|'^''^+S|^^'^^ 

+1:^'^^"^ [n, s, O + sl'^''^' + (181) 

where the first term is given by ()148p or ()149p depend- 



ing on the operator studied, and the second, fourth, and 
fifth is obtained from (dH]), as weU as P^ - P^ 

using the replacement J Jc- For the evaluation of the 

remaining contribution S^'^^'^'' |a=o from (|135p let us first 
consider the operator B = . Using 

^ L^'' P,{z,)b''+^ Pj{z,)L^'' <x (182) 

k—x.y,z 

for j = 1 and i — 0,1 (in all other cases the l.h.s. van- 
ishes) wc find 

S+i^'°^(r!,5,e,e')«i'- (183) 

We note that in (|182p the zero eigenvalue of the effec- 
tive Liouvillian appears in the resolvent left to As 
we will show in the next section the resulting term (|183p 
does not contribute to the correlation functions. For the 
evaluation of the kernel for the calculation of the suscep- 
tibility we use 



Y: i-P.(z.)Si")p,(.,)i^^^ = |f^^'^""^^'^'^'^'^^' (184) 



k—x^y,z 



0, otherwise, 
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which results in 



A=0 



min + h)- n2{n -h)- + h) + mi^ - h)) (i" - ^ i") 
I (n2in + h) + n2{n -h)- n2{i + h)-n2{i- h)) l'^ 



JrJl 



2 17 - e + i((5 - e) 



n2{^ + v + h) + n2{n -v + h)- n2i^ + v-h)- n2{n -v-h) 

~'H2{^ + v + h)~n2{^-v + h) + n2{^ + v-h) + 7^2(e -v-h)) - i l") 



(185) 



-- [n2{n + v + h) + H2{^ -V + h) + H2{n + v -h) + H2{^-v -h) 
-n2{( + v + h)- 7^2(e -v + h)- n2{( + v-h)- ■H2{i -v-h)) l'^ 

I 

Next we consider the case B = and start with the rate Fi in the second term. As can be easily shown this 
evalution of the kernel E^^^"^ from (|135p by using 



L'''P,{z,)B^^^P,{z,)L 

k—x^y.z 



2k 



±4 



O£o) r5 



, 0, otherwise. 



(186) 



sum vanishes in second order, i.e. 

S+f")(17,<S,e,0 = 0(J,3)- 
Similarly, the kernel is evaluated using 



Hence the double sum in ()135p reduces to a sum over 
i = 0, 1, where the two terms have opposite signs and 
otherwise equal each other up to the appearance of the which results in 



k—x.y,z 



'4^±' 
i r4 



(187) 



0, otherwise, 

(188) 



1 



JrJl 



{ni{n) - H2(e + h)) L% + (n2in -h)- Uiio) l 



2n-i-h + i{5 -e,') 



Hiin + V)+Hi{n-V)- H2{S. + V + h)-H2{S.-V + h))Ll 



n2{n + V-h)+ H2in -V -h) -Hi{^ + V) - Hiii ~V))Li 



(189) 



A redefinition Ac — ^ A^ yields the same result for the 
kernels with the replacement Jc J'c as can be easily 
shown in the same way as for the Liouvillian. 

Finally, let us consider the kernel in first order in the 
regime ft ^ V,ho. Starting from (|113p we obtain using 
/Ca(z) =iA/z + ... (|z| > A) 



d y.±,(l) 



1 1 Ml) 
" f7 2A^i2 



(0) 



x(e 



12 



B 



(190) 



Here the second vertex still depends on A. When i 



inte- 



grating this from max{V,ft,o} to Ac — |fi| we obtain a 
contribution cx 1/51. Hence we deduce 



v.±,(l)c ^±,(1) 



CX - for f7 > l/,/io. (191) 



A similar analysis shows that the difference between 
^^Vax{v,ho} a^^d S^''^''' is at least cx l/n. 
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VII. LONGITUDINAL SPIN-SPIN 
CORRELATION FUNCTIONS 

In this section we will use the results for the Liouvil- 
lian and the kernel to calculate the correlation functions 
([23|) and ([24|l . We first calculate the auxiliary correlation 
functions P5|) using ([M)) . We start with the parametriza- 
tion (|165p of the Liouvillian, which implies the form ((711) 
for the resolvent in C^g{Q). We further deduce from the 
previous section that the kernels T,^^ (il, iO+) admit the 
parametrizations 



Sj4r2,iO- 



3z 



+/4417)L'' + r+"(f})L", (192) 



+Y-f{n)L''+V-gf{^)L\ (193) 



where for example r^i'^'^(r2) = ftrJ^ — ^{Jr 
We stress that in contrast to the parametrization of the 
Liouvillian we have not introduced additional factors of 
i here. The stationary reduced density matrix has the 
form 



St 

Ps 



PTT 
PU 



(194) 



with + pii = 1. 
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FIG. 8: (color online) Longitudinal correlation function 
51,(0) = Ss^s^i^^) in the symmetric Kondo model (r = 1) 
for various values of the applied voltage V. 
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FIG. 9: (color online) Imaginary part of the longitudinal sus- 
ceptibility x'l (O) = Xs^ (^) in the symmetric Kondo model 
(r = 1) for various values of the applied voltage V . 



A. Longitudinal correlation functions without 
magnetic field 

The stationary reduced density matrix can be deter- 
mined using (j22p . Without magnetic field one simply 
finds = = 1/2. Furthermore, the rates (|172p and 
(|173p are given by 



Ti = r. 



TT Jb Jl v. 



(195) 



We now rewrite the resolvent using the projectors (I7ip 
and use 



Ps* 



Es-(^^)^^(^^)s+ (r!,io 

Ss^(r!)P,(r!) ^s.{n,iO+)p's' 



1 



(196) 



Tr5 
Trs 



where we have applied ([68|1 as well as (|66|) . We note in 
particular that the term rJf(rj)L'^ does not contribute 
to p96l) . This yields with jMl) and \i{VL) = -ir°(f^): 



^^J{JR + JLf5^l, (197) 



Cg^g^ (fi) 



Since 



|/mr°(r2)| 



(198) 
(199) 

(200) 



we can neglect 7mr°(ll) in (|198p . On the other hand, 
the real part of in the denominator has to be kept as 
it becomes large compared to J7 in the small-frequency 
limit. Hence we arrive at 



1 



Re F"(n) 



2 + ReV^in)"^' 



-jiJR + JL? 



ReT''{n)+in 
n^ + ReT^ifl)^ 



(201) 
(202) 



We note that the leading term of the susceptibility is 
of order J^. The correlation functions are plotted in 
Figs. [SHTOl We observe very good agreement with the 
results obtained by Fritsch and Kehrein using the flow- 
equation metho d^^i^^ . 

Let us further study the behavior of the correlation 
functions analytically. For small values of the frequency 
Ss^S^i^) has the Lorentzian form 



1 Ti 



2 f72 + F2 



Q^O 1 

2 ^ 2r^ 



1 



271 JrJlV 



(203) 



where we have used i?eF''(0) = Fi. This result for the 
small-frequency regime agrees with conclusions drawn 
from a mapping of the spin correlators to the one-particle 
Green's function of Majorana fermion o^^i^^ . On the other 
hand, in the limit of large frequencies we find 



TT (Jr_ 

4 



oc 



1 1 

-i K 



(204) 
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FIG. 10: (color online) Imaginary part of the longitudinal 
susceptibility xs^s^(^) for ^ — 100 and various values of 
the asymmetry ratio r = Jl/Jr- 
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FIG. 11: (color online) Longitudinal fluctuation-dissipation 
ratio /_L (f2) for various values of the asymmetry r and applied 
voltage V . In order to get a smooth behavior at ~ 1/ we 
have kept the arctan in the definition of Tii in this region. 
The dotted line is a guide to the eye. 



in agreement with the flow-equation method^^i^. We 
note that the J' s appearing in the correlation function 
(|20ip have their origin in the resolvent l/(f7 — L^-^ (il)), 
hence the external frequency f2 serves as a cut-off param- 
eter in Ac. This results in the logarithmic corrections at 
large frequencies. 

The susceptibility in the limiting regimes reads 



(l + r)2 



oo. 



(205) 
(206) 



The first result shows a dependence of the gradient at 
small on the asymmetry ratio r = Jl/Jr, while the 
second result indicates the revival of the fluctuation- 
dissipation theorem ([35|l for ^ y. We note that the 
derivation of (|206p relies on the fact that the coupling 
constants Jc appearing in the kernel E^. are cut-off by 
the external frequency f2, as it was discussed at the end 
of Sec. I VII Furthermore, the susceptibility x's'S'i^) 
a maximum at w Fi, where it takes the value 



^rV 



(207) 



This behavior was also deduced using the flow-equation 
methodSS.. 

In order to investigate the revival of the fluctuation- 
dissipation theorem we introduce the longitudinal 



FIG. 12: (color online) Longitudinal correlation function 
Ss'S'iyi') in the symmetric Kondo model (r = 1) for V = 
100 Tif and various values of the applied magnetic field ho. 
For ho — 115 Tk we already have V < h, which implies 
Ss'S'i^ < ft - V) = in order (see SeclVnCj. 



fluctuation-dissipation ratio^'^'^ 



(208) 



which is in equilibrium simply given by /^(r^) = 
tanh^ sgn(f^) (T ^ 0). Using our resuhs (j20T|) 
and (|202p we obtain 



(209) 



which is plotted in Fig.[nj We find /^(f} >V) = 1, i.e. 
the equilibrium result, whereas for small frequencies we 
get 



fL{^<.V) 



(l + r)2 n 



(210) 



in agreement with Refs. I2dl21l . We note that fL{^ < V) 
increases with increasing asymmetry r as the coupling of 
the voltage to the dot becomes less effective. 



B. Longitudinal correlation functions in a weak 
magnetic field {V > h) 

In the presence of an external magnetic field the sta- 
tionary reduced density matrix is given by p|| = 1 + M, 
Pll ~ ^ ~ PH ~ Pi]' ~ 0, with the magnetization in 
leading order^ (see also Refs. |1M9.,22,.24,) 



M 



IF^^(O) 1 {l + ry 



2 F«(0) 2 {1 + r-^)h + 2rV 
To evaluate the correlation function we now use 



(211) 



Trs 



M- 



r»(n) 



0, 



0, 
1, 



(212) 
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FIG. 13: (color online) Imaginary part of the longitudinal 
susceptibility Xs^S'i^) in the symmetric Kondo model (r — 
1) for V — 100 Tk and various values of the applied magnetic 
field ho. For ho = 120 Tk we have V < h, which implies 
Xs-S'{0. <h-V)=Oin order J,^ (see SeclWCj. 



where the term rg^J(fi) L'^ again does not contribute. 
From this a straightforward calculation using (|200p yields 
the correlation function up to 0{J^) 



1 Re r"(r2) + 2Af r3^(r2) 

2 n^ + Re r'^(r2)2 



(213) 



where the zero-frequency (5-peak does not appear because 
of our definition (|23p . The suppression of the correlation 
function by the finite magnetic field is shown in Fig. I12[ 
which agrees very well with similar plots obtained using 
the flow-equation mcthod'^'^. In the zero-frequency limit 
we find 



-{Jr + JlY + M ReVgfin) 



(218) 



In (|217p we have kept the terms in the second line, which 
are of 0{J^), as due to 17 /to F'f (f7) {Q ^ 0) they 
become dominant in the small- frequency limit. For larger 
frequencies these terms have to be neglected. Thus the 
static susceptibility (I32|) is given in leading order by 



r{l+rfV 



{{1 + r^)h + 2rVy 



(219) 



in agreement with the literaturei^^'^'^liS^. For larger 
frequencies the real part of the susceptibihty possesses 
logarithmic features at = h,V zt h due to the term 
ImTgf^{Q). For example, 



{n - hf + f 2 



(220) 



where the terms represented by the dots do not contain 
any logarithmic features at = /i. The imaginary part 
of the susceptibility is plotted in Fig. [T31 It has a finite 
gradient at = given by 



4r2 {I + r + r'^)h + rV 
ttJrJl {{1 + r^)h + 2rVf 



(V-h), 



(214) 

while the leading term cx 1/il in the large frequency 
regime is given by (j204p (including the logarithmic cor- 
rections in the coupling constants). Furthermore we 
observe a weak feature at O — V — h which has for 
r2 V — h ^ h the line shape 



5, 



2r2 



+ -(Q.-V + h) arctan 



- 4M(1 + r) 

n-v + }i 



(215) 



Similar features appear aX = h^V -\- h. 

For the calculation of the susceptibility we need 



Tr 



-r 



-3z 



which directly yields 

x's^s^m 



(216) 



Q^ + ReV^iQ)^ 
Jl^ + M ReTgf{Q)) i?er"(f7) 



MnimTgfin) 

, (217) 



'■V 



^((1 



)h + 2rVy JrJl 



(221) 



as well as a maximum at f2 « Fi, where it takes the value 



Xs'S-i^ 



{l + rY' 



2 {{1 + r^)h + 2rVy 



-■^Xs-S'- 



(222) 

In the large- frequency limit x's^s^i^ ^ V,h) coincides 
with the correlation function (|204|) . Furthermore, the 
imaginary part of the susceptibility has features at O = 
h,V ± h which have their origin in the function /to 7^2 
contained in Re rgf{n) and hence have a line shape sim- 
ilar to (I2T51) . 

The fluctuation-dissipation ratio fhi^) defined in 
(|208p reads in the presence of a magnetic field 



2MReVg''' 



ReT^{n)+n{jR + Jl^MH 



(223) 



which is plotted in Fig. (TH Larger values of the magnetic 
field push the system closer to its equilibrium behavior, as 
only those lead electrons in the energy interval V — h can 
couple to the dot and thus induce the non-equilibrium 
behavior. We note, however, that the equilibrium result 
/^(ri) = 1 is only reached ioi fl > V + h. Furthermore, 
increasing the asymmetry r drives the system towards 
the equilibrium situation as the coupling of the voltage 
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FIG. 14: (color online) Fluctuation-dissipation ratio /L(fi) 
for V = SO Tk and different values of the asymmetry r and 
magnetic field ho. In order to get a smooth behavior at f2 ~ 

V — h we have kept the arctan in the definition of Tii. For 

V < h the fluctuation-dissipation ratio is independently of Q. 
given by fhiP) ~ 1- The dotted line is a guide to the eye. 




FIG. 16: (color online) Derivative of the real part of the 
longitudinal susceptibility fl^ ^xs^s^(^) in the symmetric 
Kondo model (r — 1) for ho = 100 Tk and various values of 
the applied voltage V. We observe characteristic logarithmic 
features at Q = h, h i V . 
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FIG. 15: (color online) Imaginary part of the longitudinal 
susceptibility Xs»s^(^^) in the symmetric Kondo model (r = 
1) for ho = 100 Tk and various values of the applied voltage Y . 
The correlation function is given by 53252 (fi) = xs^s^(^^). 
For ^ <h — V the susceptibility vanishes in order . The 
line shape close to f2 = /i is given by (I225|l . 



to the dot becomes less effective. This effect is suppressed 
by increasing the magnetic field as overall less electrons 
couple to the dot. For small frequencies we obtain 



1 



2 (l + r + r^)h + rV V ~ h 



(224) 



C. Longitudinal correlation functions in a strong 
magnetic field {V < h) 

In the case of a strong magnetic field, V < h, the cor- 
relation functions up to quadratic order in the coupling 
are still given by ^1^, (PTT)) and (pi^ . respectively, 
where the magnetization is simply M — —1/2. One 
can easily show using /to 7^2 (f^) = fl^^l that Xs'S'^i^ 
0) — Ss's^i^ > 0), which implies the equilibrium result 
/^(ri > 0) — 1. Furthermore, the correlation function 
vanishes identically in order for fl < h — V . Phys- 
ically the Kondo spin is in its ground state \l) and the 
energy difference to the state |t) due to the external mag- 
netic field is given by h. Hence one has to apply at least 
the frequency h~V to obtain any response from the spin, 
where the energy V is provided by the applied voltage. 



This has to be contrasted with the result for the sus- 
ceptibility in the equilibrium Kondo model derived by 
Garst et al.— . They used a relation between the in- 
elastic electron scattering and the correlation function to 
show that the susceptibility in equilibrium has the small- 
frequency behavior x's^s^^^) ^ ^t^^ i-^- it is non-zero 
for < h. This linear behavior was also observed by 
Costi and Kieffer— as well as Hewson^ using a numer- 
ical renormalization group calculation. In analogy, we 
expect the non-equilibrium correlation functions to be 
nonzero for U, < h ^ V in higher order in J^. The con- 
sistent calculation of terms ~ in the real-time RG 
procedure applied here would involve, however, 5-loop 
diagrams and is hence beyond the scope of this work. 

The correlation function in the regime V < his plotted 
in Fig. [TS] We find excellent agreement with numerical 
results recently obtained by Fritsch and Kehrein using 
the flow-equation method^. In particular, we observe a 
splitting of the sharp edge at = /i due to the applied 
voltage, which leads to characteristic features at = 
h,hzt V. Using our result (|213p we can derive analytic 
expressions for the line shape close to these frequencies. 
For example, at f2 « wc find 



ttJrJlV tt 2^-/1 



4f]2 



n-h 



arctan ■ 



(225) 



The first term shows that the gradient of Xs2s^(^) ^^^^ 
become negative for Q, < h ii the applied voltage is large 
enough, i.e. V > h/2. In the vicinity of 17 = h zL V 
the correlation function shows similar kink-like behavior 
(|225p . The physical origin of these kinks lies in the fact 
that at each of the energies ^ = h, h zt V a new process 
sets in, which involves a spin-flip on the dot costing the 
Zeeman energy h as well as the virtual hopping of an 
electron on and off the dot gaining or costing the energy 
—V, 0, or V, respectively. The real part x's^sA^) of 
the susceptibility shows logarihmic features (|220p at = 
h, hztV as is shown in Fig. [111 We stress that the splitting 
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FIG. 17: (color online) Imaginary part of the transverse sus- 
ceptibility Xt(^) = Xs-s+ (^) symmetric Kondo model 
(r = 1) for ho = 100 Tk and various values of the applied volt- 
age V. 




FIG. 18: (color online) Imaginary part of the transverse sus- 
ceptibility x's-s+(^^) for V = 200 Tk, /lo = 100 Tk, and var- 
ious values of the asymmetry ratio r — Jl/Jr- The result is 
invariant under r 1/r. 



of the sharp edge at 
effect. 



h is a true non-equihbrium 



VIII. TRANSVERSE CORRELATION 
FUNCTIONS 



FinaUy let us discuss the transverse correlation func- 
tions in the presence of a magnetic field. We note that 
by virtue of we can restrict ourselves to the S~S'^- 
correlations. The corresponding kernels up to second or- 
der in Jc were calculated in (fTig)) . ([115)) . 
(fTS5)) -- (fT35)) as well as We will first discuss the 

susceptibility and present the results for the correlation 
function afterwards. 

In order to derive the susceptibility we start with the 
parametrization 



E5+(r!,io+) = r-2L2 



r-3it + rs;'(f^)it> (226) 



where for example F 



-2i - itrJc 



-2i - + Jl) + + Jl) ■ Note that we already 

indicated that the explicit frequency dependence in or- 
der appears in r^|^ exclusively. (There is of course an 
implicit frequency dependence of and TZ"^ through 



Jf,.) Now using 



Trc 



we obtain 



I]s-(r!)P,(17)I],+ (r!,iO+)p^* 



i {t--^ + {T--^ + 2r-;'(17))Af) S^+ (227) 



- s+ 



(r-;V2r^;^(r!))A/ 



VL - h{VL)+iT2{n) 



(228) 



where we have introduced the short-hand notation 
r2(0) = r°(f])-hr^(f^). For /lo = Owe find xs-s+(f^) = 
2x5^Si(ri). The transverse susceptibility has a peak at 
the solution of 



- i?e /i(r2) - /to Fa (r?) = 0, 

which is up to first order solved by 
1 



1 



-AJr + Jl) ]ho = h. 



(229) 



(230) 



In a finite magnetic field the spin on the dot will be in 
its ground state ||). The energy difference to the excited 
state It) is given by h, leading to an enhanced response of 
the system at this frequency. At the peak the imaginary 
part of the susceptibility takes the value 



1 2 



(231) 



as is shown in Figs. [T7I and [T51 The peak is suppressed by 
increasing the voltage, since this reduces the probability 
for the Kondo spin to be in its ground state. On the other 
hand, for a fixed value of the voltage the peak increases 
with increasing asymmetry ratio as the coupling of the 
voltage to the dot becomes less effective. The width of 
the peak is up to order given by 



ReT2{h)- Imh{h) = f 2 



with the limiting cases 



h^V: ttJrJlV. 



(232) 

(233) 
(234) 



In the equilibrium limit, V = 0, this corresponds to the 
result obtained in Ref. [65|. The real part of the trans- 
verse susceptibility possesses logarithmic features similar 
to (|^^ a.tn = h,h±V. 

Using a pseudo-fermion representation of the Kondo 
spin together with non-equilibrium perturbation theory 
Paaske et al.— previously obtained the transverse suscep- 
tibility. In order to compare these results to (|228p we first 
make the approximations h{n) ho and F2(ri) Fa- In 
the limit ho ^ we then obtain using F^4 = ^(Jr+Jl)'^ 



(235) 
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In the regime Hq ^ f 2 we have M = 0{J^), thus we can 
neglect F^l as well as r~4 in the numerator in 
which results in 



Xs-s+i^) 



2M 



Q - ho + iT2 



(236) 



These approximations agree with the results obtained in 
Ref. m (we have to replace Hq ^ —B due to a different 
sign in the definition of the bare dot Hamiltonian Hs). 
Furthermore, in the regime fl^ha < V we can use 



Af 

's+ + (^s;" + "^^s^'i^)) M -2— fs - 2iM, (237) 



where we have replaced ImTii{Q) in tbc real 

part of r~|^, to obtain 



Xs-s+i^) 



2M ho - iFs 
^0 fl- hQ + if 2 



(238) 



This confirms a conjecture by Paaske et ali^. We would 
like to stress, however, that our result (|228p goes beyond 
the approximation (|238p . 

In analogy to the susceptibility one finds for the cor- 
relation function 



%M{jR + jLf{n-ho) 



{n-Reh{n)-Im T2{n))^ + {Re r2(o)- 



-Imh{n)f' 

(239) 



where we have neglected all terms of order in the 
numerator. This allows the calculation of the transverse 
fluctuation-dissipation ratio 



'>s-s+ 



(240) 



which is plotted in Fig. 1191 For negative frequencies 
H. < —V the fluctuation-dissipation ratio takes the value 
fxi^) = —1, whereas for frequencies H. > V we flnd 
fri^) — Ij thus recovering the equilibrium situation 
in these limits. As for the longitudinal fluctuation- 
dissipation ration we observe that increasing the mag- 
netic field or the asymmetry ratio r drives the system 
towards the equilibrium situation. 



IX. CONCLUSIONS 

In this article we have generalized the real-time renor- 
malization group method in frequency space to allow the 
calculation of dynamical correlation functions of arbi- 
trary dot operators in systems describing spin and/or 
orbital fluctuations. We applied this to the two-lead 
Kondo model in a magnetic fleld, where we calcualted 
the longitudinal and transverse spin-spin correlation and 
response functions up to second order in the exchange 
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FIG. 19: (color online) Transverse fluctuation-dissipation ra- 
tio fri^) for V = 80 Tk and various values of the asymmetry 
ration r and the applied magnetic field ho. The plot shows 
good agreement with similar results obtained in Ref. |6l|. The 
dotted line is a guide to the eye. 



coupling. We wish to stress that within this formal- 
ism the Kondo spin is directly represented by matrices 
in Liouville space, hence there is no need to apply a 
pscudo-fcrmion representation. Speciflcally, we derived 
the two-loop RG equations for the dot operators and 
solved them analytically up to order in the weak- 
coupling regime. Here Jc denotes the effective coupling 
at the energy scale Ac = max{y, ho} which has to satisfy 
Ac ^ Tk- Our results show several features attributed 
to the non-equilibrium situation, e.g. the splitting of the 
edge at i7 = ft. of the longitudinal correlation function 
in a strong magnetic field or the suppression of the peak 
in the transverse susceptibility by a finite applied volt- 
age. Furthermore, we find very good agreement with 
results for the longitudinal correlation function recently 
obtained by Fritsch and Kehrein using the flow-equation 
metho d^^'^° . A particular advantage of our approach is 
the possibility to obtain analytic expressions for all cor- 
relation functions in the weak coupling limit. 

We have calculated the spin-spin correlation functions 
for the nonequilibrium Kondo model in the weak cou- 
pling regime Ac S> Tk- The regime of strong coupling, 
Ac < Tk, is still an open problem. In this case, the 
exchange couphngs Jc become of order 0(1) and a con- 
trolled truncation of the RG equations is no longer possi- 
ble. Within the present RTRG-FS method it was shown 
in Ref. [l| that the relaxation/dephasing rates saturate to 
the Kondo temperature in the strong coupling regime. 
As an effect the coupling constants do not diverge as 
in poor man scaling methods but remain flnite. How- 
ever, the numerical solution of the RG equations in low- 
est order showed an instability against an exponcntion- 
ally small change in the initial condition for the relax- 
ation/dephasing rates. Although it was possible to flnd 
excellent agreement for the temperature dependence of 
the linear conductance with NRG calculations, it was 
necessary to fine tune the initial condition for the rates. 
Therefore, up to now, it is not yet clear whether a con- 
trolled solution of the strong coupling regime is possible 
by using RTRG-FS. 

The nonequilibrium Kondo model describes the spin 
fluctuation (or Coulomb blockade) regime of the more 
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general nonequilibrium Anderson impurity model (for a 
systematic derivation of the Kendo model from the An- 
derson model using a Schrieffer- Wolff transformation, see 
e.g. Ref. 37). In this model the single-particle spec- 
tral function is of most interest, which has recently been 
studied within NRG in a scattering wave basi a^^i"^" . The 
present RTRG-FS method can also be applied to this 
model but, as explained in detail in Ref. [l|, in the charge 
fluctuation regime it is not yet clear whether a well- 
defined weak coupling regime exists at zero temperature. 
At resonance (i.e. when the renormalized single-particle 
level is identical to one of the chemical potentials of the 
leads) there is no energy scale except the broadening of 
the level itself and the expansion parameter is of order 
0(1). Nevertheless the results obtained in Ref. [s^ (using 
a previous version of the real-time RG method) after an a 
priori uncontrolled truncation of the RG equations were 
in excellent agreement with the Bethe Ansatz solution^ 
for the equilibrium occupation of the local level. These 
and related topics are of high interest and will be the 
subject of forthcoming research. 
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APPENDIX A: DERIVATION OF JSO]) AND (f3T|) 

Obviously the time-dependent correlation function can 
be written as 



SABit) - i([^(t)H,S(0)H] 



Now, applying the time-translational invariance 



Aih)iiB{t2)R) = (A{h+t)nB{t2+t)u) (A2) 

' St \ I St 



for all t, ii, t2 fixed and finite as well as the relation 



[v4(0)H,i?(t)H]±) = ([v4(t)H,i?(0)H] 



(A3) 



which can be verified by a straightforward calculation 
using i? = At , we obtain pO]) . In the same way (jA2[) can 
be used to derive the relation xab^) — iC'As(^)- 



APPENDIX B: DERIVATION OF B^^^^, 



(2) 

In this appendix we will calculate the vertex Oj^ -^y. 
Starting from ((96|) we first take the zero-temperature 
limit. Furthermore, we can expand the resolvents up to 
0(J) as38 



n(^,c^) 



E^iuj- Lf - l';^ -{E + itj)Z(i) 



= 1- 



2 } E + \u:-Lf-Lf 



1-^1 (Bl) 



with 



_ r(l) 1 



^ 2 



^(1) ^(0) 



(B2) 



{A)st{B)sf (Al) 



Using this we obtain together with the expansion B± 



d - ( 1 

— 6±,ii'(f7,(5,^,C';wi,t^i') = 1012 1 TT M — ^ 

V 2 y + 15 + iA + iwi - Lf - L^s' 



1 



^2 + + iA + ic^i - Lf - 



1 - 



G2i,-(1^1'), (B3) 



where we have omitted the arguments of the vertices G for simplicity. Using 



L 



(0) 



L 



(0) 



riO) , f(l) 



^(l)^g(O) 



^(l)^g(O) 



(B4) 
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with k = zb/ig T tr J for B = and k — for B — (see Sec. IVI[) we obtain after a partial fraction expansion 



1 



r2 - C - K + i((5 - 



Gl2 1 



B. 



(0) 



(1) 



1 r 
2 



57i2 — K + i5 + iA + ibJi — L 



(0) 



$12 + if + iA + icji - L^°^ - L^'^ 



G21. + (1 ^ 1'). (B5) 



We see that no terms cx 1/A can occur on the r.h.s. and, hence, in order to determine the RG equation for B^\y we 

have to replace the vertices G by the leading order ones G*^^^ and omit all terms containing /B^-*, Z^^'^ and l'^q^ ■ This 
yields using k~K + 0{J) as well as j^^^ = ^ + In ^^^^ ■ 



A ' rfA 
1 



SI - C - K + i((5 - f ) 



dA A 



B. 



(0) 



A + ^i -i(Ci2+i$' -^ 
A 



(0)^ 

S ' 



G^i^, + (1 ^ 1'). 



(B6) 



ii , 11 il li U , il 

2 3 1 r 3 2 2 3 1 1' S±3 2 



1 2 g± 1'3 3 2 1 2 r3 3 2 



Irnl ll li 1 1 . ll 



2 3 3 1 2 1' 2 3 3 1 S± 2 1' 
FIG. 20: (color online) Two-loop diagrams for B±^iii 



1 ^ ^ 1^^ 

12 2 1' 12 2 1' 

B±,12 S±,2l' 

FIG. 21: (color online) One-loop RG diagrams for 
which contain B± n' itself. 



Using the RG equation for the leading-order vertex (|101[) 
we can integrate (|B6p up to higher-order corrections and 
obtain the solution (jllSp , which in particular satisfies the 
initial condition 



4'n'(f^>'5,C,C';a;i,c.iO 







A=Ao 



(B7) 



as given by ([75]) from the discrete RG step. 

The two-loop diagrams for the vertex B±^iii are shown 
in Fig. [201 As each diagram contains three vertices G (x J 
and we are interested in the vertex up to second order, 
we have to extract the terms ~ 1/A. This is done by 



1 



expanding the resolvents in lowest order 

Il{E23,u; + A + L03) « --r—, (B8) 
U{Ei2,uj + A + LJi) « --, (B9) 



n(-Ell'23, LJ + A + LUl + LUy 



A + UJ3' 

i 

"A' 

i 



, (BIO) 



A + LU3 

where E (oj) stands for either (6) or ^ (^'), respectively. 
Using this we immediately see that the first two diagrams 
are proportional to 



du}3 



the fourth and fifth to 



1 

A2' 



1 



dLLi3 



A Jo {A + UJ3) 
and the third and sixth to 



2 °^ ^2' 



1 

A2 



duJ3 1 

oc 



A -I- W3 A2 



(Bll) 



(B12) 



(B13) 



Thus the two-loop diagrams behave as ~ J^/A^ and 
hence do not contribute to the renormalization of the 
second-order vertex B^^\y. 

Finally we have to study the one-loop diagrams which 
contain B±^iii itself. These diagrams are shown in 
Fig. [21] As we can easily see from (jllSI) . the leading- 
order result for B±^ii' behaves for large A as 



B^^\,,{n,5,^,e;^i,^2) ^ jG' - ^. 



(B14) 



Now the one-loop diagrams contain an additional vertex 
G as well as a resolvent 11 ^ 1/A. Hence we deduce that 
the diagrams yield terms proportional to / A^ which do 
not contribute to the renormalization of the second-order 
vertex 
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APPENDIX C: DERIVATION OF S^'^^^ 



In this appendix we will derive the RG equation (|119p 
for the second-order kernel T,^ {Q, S, ^, ^'). For this we 



have to evaluate all terms on the r.h.s. of ([97]) . We start 
with the first line. Using (jBip and expanding G in powers 
of J yields 



( 1 



2 ;j7i2+i<5 + iA + ic.2-4°'-4^^ ^ 



1 r 



^(l)^g(O) 

Z(i) 



- ^dcu^ (iG^^^+Gif (f},<5;A,L.2)) 
Jo 



■ ^12 + if + iA + icj2 - Lf ~ il!^ 



(1) 



2 /^2i 



g 

+ i(5 + iA + icj2 - is ^ 



(0) i 

62 + i^' + iA + icj2 - 



(CI) 



G (C2) 



aLL>2 'j-12 TTrT "± 77TT 

r!i2 + 1(5 + iA + icj2 - ' 62 + if + iA + ic^2 - Lf 



X (iG^j"^) + (5^/^(62, e' + A + C.2; -CU2, -A)), 



(C3) 



where we have already neglected the terms Z^^'^ and L^^-* in the second and third line as they lead only to higher-order 

corrections. Using the commutators (jB4p the first line (jCip can be treated similarly to S^'*^^^ derived in Sec. IV Cl 

In the term ex yB^"* we only keep the term containing ICAiz) (the term cx 1/A was already used to calculate S^'*'^^), 
while in the other two terms we have to extract the term cx 1/A. Thus we arrive at 



G 



(1) 



/CA(r!i2 + i^ - 4"^)sf - /Ca(62 + - 40 



(0)^ 



G 



(1) 



^ 2A ^2 =^ 21 + :JX<^12 
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^(1)^^(0) 



G 



(1) 



(C4) 



Using the same steps for the terms containing G^^"^-' in (jC2p and (|C3[) one finds 



1 

2A 



■^(2ai) ^,(0) ^(1) ^(1) K,(0) ^(2ai) 
•-^12 ''i <-^21 +'-'12 ^± '^21 



(C5) 



For the terms of (IC2| and jCS]) containing G^'^^^ we use ([TU5)) in the form 



Gi2 [^,0;A,uj2)^Gy In 



^32 ~ ^23 



(1), A + c^2 - i(»23 + i'^ - 4°0 



G^il (C6) 



A "32 "^3 ^ 

G^f (C12, + A+C.2; -C.2, - A) . G« In '^-^(^"+^^^-^^°') G(^,) - G^^ In ^ + ' ^^^^^ + ' ^^"'^ G« . (C7) 



A 



When inserted into (jC2p and (|C3p the first terms do not depend on the integration variable lo2- The remaining integral 
can be done as usual by a partial fraction expansion. This yields (k = ±ft,o foi' S = and /t = for B = S^) 



^ rW 2A-i(»i3+i<5-40 Ml) 



/CA(f^l2 +16- Lf) - /CA(a2 + ie' - i^'O 
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(1) 
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/CA(r!i2 -f 15 - 4^6^ - /CA(ei2 + i^' - L's') 
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-(0)n 



(1) 2A- 1(^13 + 1^-40 ^(1) 



If we now expand /Ca and the logarithm for large A, 
/Ca(z) = ln2 + iz/2A and In = ln2 - iz/2A, and 



keep only the terms proportional to ^ we arrive at 



In 2 
2A 



G 



'^23 ^± 



CT23 



g: 
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31 ■ 



(C8) 



29 



In contrast, the second terms of (|C6|) and (|C7|) do de- 
pend on the integration variable W2- The evahiations is, 
however, straightforward. We use a partial fraction ex- 
pansion for the resolvents left and right to the vertex B^'' 
as well as 



In 



A + UJ2 ~ iz 



In 



A + UJ2 



IZ 



A A A + W2 

This leads to integrals of the form 



duJ2 —, : — In , 

z + 1A + 1UJ2 A 



A + UJ2 1 - In 2 



2A 



(C9) 



z, (CIO) 



dLU2 7 : 

Z + 1A + ILU2 A + LU2 



-2A' 



where we are only interested in the ~ 1/A terms. Now 



using the asymmetry 

^(1 



-G^i we find 



2A 



Ct23 i5± 



G 
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(C12) 



which has to be combined with (jCSp for the full result 
from the terms containing 

The second and third line of ([57)1 containing the vertex 

(2) 

Bj_\y can be treated using the same steps as were used 

to evaluate the G'^^''^ -dependent parts of ((C2)) and l(C3)) . 
The result reads after some tedious but straightforward 
algebra 



1 + ln2 



/^(i) /^(.J-; RV"; I K!\") r<^ 
^12 "'23 ± ± 2 
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23 



(C13) 



Finally, to extract the leading term ~ J^/A oi the 
fourth line of ((97|) one can simply replace the resolvents 

U{z,z' + A + uJi) by 1/{A + uji). This yields 



In 2 
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G\ 
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•^23 "± 



g: 



(1) 

31 ■ 



(CM) 



Hence, the result for the RG equation of the kernel 
E|'^^^(f7,<5,e,C') is obtained by summing jCl, jCS]), 
del), dCllI, jClll and den, using ]Ca{z) = ^Fa(z) 
with (fnel) in (|C4| as well as 



X^12 



^(1) 0(0) , ny^j ^yi-) ^y-) _ /^y-) n\^) /^y-. 
•-^23 "± ~^ "± ^23 ^31 ~ ^ ^12 "± *-^21 



3(0) ^(1) 



This yields pl^)) . 



(1) r(0) ^(1) 
^21 ■ 

(C15) 



APPENDIX D: ALGEBRA IN LIOUVILLE SPACE 

Consider an operator A acting on the dot Hilbert space 
having matrix elements Aat with respect to the basis 
(It) , li)}. UK denotes the superoperator acting on dot 
operators via O. — [A, .]_^ then for an arbitrary dot op- 
erator B we have 

{OB)ab — Oab,cdBcd, Oab,cd — AacSbd^Sac Adb- (Dl) 



Furthermore, we represent superoperators in the matrix 
representation 



O = {Oa, 



( Ott^tt 



b.cd) 



OTT,ii 
Oii.Xi 



C*Ti,TT 



OiT,Ti 



On, XT 



(D2) 

If O = PQ is the product of two superoperators, then 
Oab,cd — Pab,ef Qef,cd and the matrix (|D2[) of O is simply 
given by the matrix product of the matrices of P and Q. 

A basis for the operators in the Liouville space of the 
Kondo dot can be built up by the spin superoperators 
lA and L~ defined in p38p . An explicit representation 
in the basis (|D2p is provided by 
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Furthermore we define the operators 
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L" 




(D4) 
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= il + 2L+^i-^ 


(D5) 
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= \{L+ -L-)-iL+ xL-, 


(D7) 
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= ^(i^ -JT) +^L.+ X IT, 


(D9) 


as well as 






LI = 


L^^±iL''y, a = +,-,1,2,3, 


(DIO) 


L% = 


L\±{L+L-'-+L+'- L^), 


(Dll) 


Ll = 


lItIl+l-^ + l+^'Li). 


(D12) 


We note that L 


= -2L2^ as well as L%+L% = 


2Ll. 
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In the spin sector we will use frequently 



<i^3^i^3<T2 = '^Izoya^a^ = 5oi<y2 (no sum over a) 



(D13) 



as well as 

'^li<yi<^l^a2 = -<,<T2'^^i<T3 '^'^^abc (tI,„^ (for a ^ h). 

(D14) 

The Liouville operators satisfy (the sums are over i,j = 
X, y, z while the index p takes the values p = ±) 



id 

i 
i 



r2k 



2i 
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^^2. J^(i) 
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Jp-^trJL''. (D19) 



where we have used zq — 0, zi = 0{J^) and z± = ±/io + 
0{Jc). Furthermore, we can expand the zero-eigenvalue 
projector as 



^o(^i)-^o(o)-i^|n.) 

= Poizo) + 0{Jc), 



(E2) 



where -^Viz) 



z=0 



-i(Ji? + Jl) + ■■■■ This directly 



yields (|133p . To prove the final statement, k 
0{Jc) for all pairs («, j) for which 



(E3) 



APPENDIX E: PROOF OF (fT32l) AND (fTM]) FOR 
THE ISOTROPIC KONDO MODEL 

The results for the Liouvillian presented in Sec. IVI Bl 
allow us to obtain 

Lf ^ ho L''^ hoiP+ - ^ J2 z,P,{zi) + 0{J,), 

i=0,l,± 

(El) 



is non- vanishing, we note that for B — these pairs 
are given by (see (fT82l) and (fTM]) ) = (0,1), (1,1) 

for S'^^' and («,j) = (±,±) for b'"^\ respectively. On the 
other hand we have k = in both cases, which implies 
Zi — Zj = K+0{J^) for B = . The same analysis can be 
performed for B — where the relevant pairs are given 

by (*, j) = (0, T), (1, T) for and (z, j) = (1, t), (±, 1) 
for B^^\ 
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